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ABSTRACT 


The  photorefractive  effect  is  a  nonlinear  optical  phenomena  that  can  produce 
phase  conjugate  (complex  conjugate)  reflection  of  visible  light  in  a  crystalline  solid. 
Very  little  is  known  about  the  microscopic  origin  of  photorefraction;  consequently  this 
thesis  investigates  the  following  microscopic  parameters  of  barium  titanate  (BaTiO^  ): 
the  quantum  efficiency,  charge  transport  length  and  the  carrier  generation  rate  both 
theoretically  and  experimentally.  BaTiO-j  wras  chosen  because  it  is  a  highly  nonlinear 
crystal  with  large  electrooptic  coefficients  and  therefore,  is  a  very  efficient 
photorefractive  material. 
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I.  INTRODUCTION 


Photorefraction  is  the  change  in  the  index  of  refraction  of  an  optical  material 
under  intense  laser  illumination.  Photorefraction  was  originally  discovered  in  nonlinear 
optical  investigations  of  LiNbO-j  [Ref.  1].  At  that  time,  it  was  interpreted  as  optical 
damage  in  the  crystal  due  to  undesirable  scattering  and  decollimation  of  laser  light. 
Later,  Kukhtarev  et  al.  [Ref.  2]  explained  it  as  a  three  step  process:  1)  Under  laser 
illumination,  optical  excitation  of  electrons  and  holes  occurs.  Electrons  move  from 
donor  sites  in  the  crystal  to  the  conduction  band  and  holes  move  from  acceptor  sites  to 
the  valence  band.  2)  The  photo-excited  carriers  migrate  from  the  illuminated  regions  to 
the  dark  regions  under  the  influence  of  diffusion,  drift  in  an  external  electric  field,  and 
the  bulk  photovoltaic  effect  (BPVE).  The  charge  transport  mechanisms  separate 
charge  pairs  from  each  other  and  this  produces  a  space-charge  field,  inside  the  crystal. 
3)  This  induced  space-charge  electric  field  alters  the  index  of  refraction  of  the  crystal 
via  an  electrooptic  (Pockel)  effect. 

If  two  or  more  beams  intersect  to  form  an  interference  pattern  within  a  crystal, 
the  standing  wave  pattern  can  produce  a  three  dimensional  volume  hologram.  This 
grating  can  then  produce  a  phase  conjugate  reflection  (complex  conjugate  of  the  E-M 
field)  of  the  incident  beam  that  is  important  in  phase  conjugation  studies.  In  addition, 
the  grating  is  responsible  for  a  transfer  of  energy  between  the  two  intersecting  beams. 
This  process  is  called  two  beam  coupling. 

Photorefraction  depends  on  the  following  parameters:  a)  the  amount  of  light 
absorbed  by  the  crystal,  b)  the  density  and  type  of  impurities  in  the  crystal,  c)  the  type 
of  carriers  and  the  charge  transport  mechanisms,  d)  and  the  magnitude  of  the 
electrooptic  coefficients. 

BaTiO-j  was  chosen  for  investigation  since  it  is  a  ferroelectric  crystal  that  has 
large  nonlinear  electrooptic  coefficients.  Consequently,  BaTiO-j  has  a  large  diffraction 
efficiency  that  makes  it  ideal  for  phase  conjugation  and  photorefraction  applications. 

Photorefractive  materials  such  as,  BaTiO-}  ,  LiN'bO-j  ,  LiTaO-j  and  bismuth 
silicon  oxide  are  finding  many  applications.  They  can  be  used  as  optical  phase 
conjugation  elements  and  as  volume  holographic  optical  elements.  Optical  phase 
conjugation  is  used  to  correct  static  and  dynamic  distortions  in  optical  systems.  Since 
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the  phase  conjugate  beam  is  a  complex  conjugate  of  the  incident  beam  wave,  it  is  in 
effect  a  time  reversed  replica  of  the  incident  beam.  Consequently,  phase  distortions 
imposed  by  propagation  through  a  perturbing  media  such  as  the  atmosphere  will  be 
undone  when  the  phase  conjugate  wave  returns  through  the  medium  along  the  same 
optical  path.  In  addition,  phase  conjugation  is  used  to  improve  laser  resonators,  in 
atmospheric  communication  networks  and  in  image  amplifying  systems.  Since  the 
interference  pattern  represents  a  volume  hologram,  photorefractive  materials  are  used 
as  volume  holographic  storage  media  in  optical  data  processing  systems. 

In  order  to  understand  the  source  of  the  photorefractive  effect  in  BaTiO^  ,  the 
quantum  efficiency  and  the  number  of  photo-excited  carriers  produced  per  absorbed 
photon,  were  needed.  Consequently,  this  thesis  attempts  to  determine  the  quantum 
efficiency,  the  charge  transport  length  and  carrier  generation  rate  in  BaTiO^  .  The 
quantum  efficiency  was  unknown  at  the  start  of  this  thesis.  To  determine  the  quantum 
efficiency,  the  photorefractive  absorption  coefficient  was  needed,  in  addition  to  the 
pyroelectric  current,  the  bulk  photovoltaic  current  and  the  phase  shift  angle  between 
the  hologram  grating  and  intensity  pattern.  Chapter  II  contains  the  main  aspects  of 
the  theories  for  spontaneous  polarization,  the  ferroelectric  properties,  the  pyroelectric 
effect,  the  bulk  photovoltaic  effect,  the  semiconductor  properties,  the  charge  transport 
mechanisms,  the  electrooptical  and  photorefractive  effect  in  BaTiO-j  .  In  addition,  the 
theoretical  section  contains  an  expression  for  the  quantum  efficiency  in  terms  of 
experimentally  measurable  parameters.  Chapter  III  presents  the  experimental 
observations. This  includes  spectrophotometer  absorption  measurements  of  BaTiO^  for 
different  wavelengths  and  orientations  of  polarized  light.  In  addition  the  pyroelectric 
effect,  the  change  in  the  spontaneous  polarization,  versus  temperature  was  used  to 
measure,  indirectly  the  absorption  coefficient.  The  quantum  efficiency  was  ultimately 
determined  from  the  pyroelectric  absorption  coefficients  and  the  phase  shift  angle 
determined  from  the  beam  coupling  measurements.  Chapter  IV  summarizes  the  results 
of  this  thesis. 
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II.  THEORY  OF  PHOTOREFRACTION 


A.  THE  STRUCTURE  OF  THE  BATIO3 


Figure  2.1  Unit  cell  of  BaTiOy 

?  4- 

Figure  2.1  shows  that  a  unit  cell  of  BaTiO^  contains  a  Ba  ion  on  each  corner 
of  the  cube,  an  O^*  ion  at  the  center  of  each  face  and  a  Ti^+  ion  at  the  center  of  the 
cube.  BaTiOj  has  four  different  unit  cell  structures  for  four  phases  [Ref.  3].  It  has  a 
cubic  symmetry  above  the  Curie  temperature  118°C;  and  it  is  nonpolar, 
centrosymmetric  and  nonpiezoelectric  with  point  group  m3m.  In  the  cubic  phase, 
BaTiOj  is  not  ferroelectric,  but  has  pyroelectric  properties  and  it  is  paraelectric.  The 
crystal  has  tetragonal  symmetry  between  118°C  and  5°C  with  point  group  4mm.  This 
phase  is  polar  and  the  four  fold  axis  is  the  polar  direction  known  as  the  c-axis.  In  this 
phase,  the  crystal  has  ferroelectric,  pyroelectric  and  semiconductor  properties,  which 
permits  photorefraction.  A  tetragonal  unit  cell  results  from  distortion  of  the  nonpolar 
cubic  cell;  one  of  the  cube  edges  stretches  to  become  the  tetragonal  c-axis,  the  other 
two  edges  compress  to  become  the  tetragonal  a-axis.  At  5°C,  another  phase  transition 
occurs  and  BaTiO^  becomes  orthorhombic  with  point  group  mm.  In  this  phase,  it  is 
still  ferroelectric  but  the  direction  of  spontaneous  polarization  is  no  longer  the  original 
cubic  <  100 >  direction  as  it  was  for  the  cubic  and  tetragonal  phases.  Instead,  it  is 
parallel  to  one  of  the  original  cubic  <  1 10  >  directions.  The  orthorhombic  unit  cell 
arises  from  distortion  of  the  cubic  cell;  one  of  the  face  diagonals  is  elongated  to 


11 


become  the  orthorhombic  a-axis,  the  two  other  face  diagonals  are  compressed  to 
become  the  orthorhombic  b-axis.  The  last  transition  is  at  -90°C  and  the  unit  cell 
becomes  rhombohedral  with  point  group  3m.  In  this  phase,  the  polar  axis  is  parallel  to 
one  of  the  original  cubic  <  1 1 1  >  directions.  The  rhombohedral  unit  cell  results  from 
the  elongation  of  the  original  cubic  cell  along  one  of  the  body  diagonals.  This  section 
can  be  summarized  in  Figure  2.2  . 


Figure  2.2  Unit  cells  of  the  four  phases  of  BaTi03  [Ref.2]. 

(a)  The  cubic  phase  is  stable  above  11S°C,  (b)  Tetragonal  is  stable  between  118°C 
and  5°C,  (c)  Orthorhombic  is  stable  between  5’C  and  -90°C,  (d)  Rhombohedral  is 
stable  below  -90°C. 

The  dashed  lines  in  Figure  2.2  (b),  (c)  and  (d)  shows  the  original  cubic  cell;  Ps  is  the 
spontaneous  polarization  vector. 
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B.  DIELECTRIC  PROPERTIES  OF  BATIO3 
I.  The  spontaneous  polarization 

In  general  there  are  three  types  of  polarization.  1)  Dipolar  polarization  arises 
from  molecules  which  have  dipole  moments  that  can  change  their  orientation  when  the 
field  is  applied.  This  type  of  polarization  is  important  in  the  LTIF  and  microwave 
region  of  the  spectrum.  2)  Ionic  polarization  results  from  the  displacements  of  the  ions 
with  respect  to  each  other  and  is  important  in  the  infrared  region  of  the  spectrum.  3) 
Electronic  polarization,  which  arises  from  the  distortion  of  the  electron  shell  relative  to 
a  nucleus,  is  dominant  in  the  optical  part  of  the  spectrum. 

In  an  anisotropic  crystal,  i.e.  the  tetragonal,  orthorhombic,  and  rhombohedral 
phases  of  BaTiO^  ,  the  polarization  is  described  by, 

?.  =  (£..-!)£„  E..  (eqn  2.1) 


In  the  last  equation,  £~  is  the  dielectric  tensor,  and 

=  ,  (eqn  2.2) 

Xjj  is  the  susceptibility  tensor.  In  the  formation  of  the  spontaneous  polarization  in 
BaTiO^  ,  some  important  concepts  are  involved  such  as,  "Lorentz  correction"  and 
"Polarization  catastrophe".  These  are  the  main  factors  that  determine  the 
ferroelectricity  of  the  crystal  (see  for  example  [Ref.  4])  In  this  presentation,  cubic 
symmetry  was  considered.  Following  Kittel  [Ref.  4],  the  Lorentz  correction  correction 
is:  the  polarizing  field  acting  on  the  crystal  has  to  be  different  from  the  local  field 

Et  =  E  +  (l/3e0)Ps .  (eqn  2.3) 

In  this  equation,  E (  is  the  local  field,  E  is  the  external  field,  P$  is  the  spontaneous 
polarization  and  £q  is  the  dielectric  permeability  of  space.  The  relation  between 
polarizability  and  the  dipole  moment  is  given  by 

P=bFj,  (eqn  2.4) 


13 


where  P  is  the  dipole  moment  of  an  atomic  or  ionic  species  (assumed  to  be  spherically 
symmetric)  and  b  is  the  polarizability. 

The  polarizability  is  an  microscopic  parameter.  The  dielectric  constant 
depends  on  the  crystal  symmetry’  which  becomes  a  tensor  for  a  non  cubic  cell.  The 
spontaneous  polarization  is  the  sum  of  all  the  dipole  moments  per  unit  volume, 

Ps  =  SNjPj  =  I  bj  E,  (j).  (eqn  2.5) 

In  the  last  relation,  Iv  is  the  concentration  and  bj  is  is  the  polarizability  of  atoms  of 
type  (j)  and  Ej  (j)  is  the  local  field  at  the  atoms  sites  (j).  The  next  step  is  to  find  a 
relation  between  the  dielectric  constant  and  polarizabilities.  Substituting  the  Equation 
2.3  into  2.5  gives  for  cubic  symmetry', 

*.-(!  bj){£  +  (fs/3  e0»  •  (eqn  2'6) 


After  some  algebra,  Equation  2.6  becomes 

(£  -1)  (£  +2)}  =  (l/3E0)£N.b.,  (eqn  2.7) 

o 

the  Clausius-Mossotti  relation.  At  optical  frequencies  n*-  =  £  and  Equation  2.7  then 
becomes 

(n2  -l)/(n2  +2)- (1/3  £0)  £  N  .  b.  (electronic)  ,  (eqn  2.8) 

where  n  is  the  index  of  refraction.  Rearranging  Equation  2.7,  the  dielectric  constant 
takes  the  form 

£  =  {1  +(2/3  £0)  £  N.  b.  }/  { l-(  1/3  £q)  y  N.  b.)  .  (eqn  2.9) 

In  the  last  relation,  bj  is  the  electronic  plus  ionic  polarizability  of  (i)  type  ion  and  N;  is 
the  concentration.  The  dielectric  constant  becomes  infinite  when  the  term  in  the 
denominator  vanishes  and  leads  to  the  "Polarization  catastrophe".  This  means  that 
BaTiO^  has  a  finite  polarization,  even  when  there  is  no  applied  field.  Indeed,  the 
catastrophe  occurs  at  the  Curie  point  where  the  dielectric  constant  as  a  function  of 
temperature  has  a  peak,  as  is  seen  from  Figure  2.5. 
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The  polarization  catastrophe  occurs  when 


(1/  3  e0)Y  N.  b.  =  1  .  (eqn  2. 10) 

In  the  original  theory  of  Devonshire  [Ref.  5],  the  ionic  species  were  assumed  to  be 
spherically  symmetric.  Slater  [Ref.  6]  showed  that  Equation  2.10  can  not  be  satisfied 
under  the  condition  of  spherical  symmetry,  so  he  introduced  asymmetry  of  the  ionic 
species. 

The  latest  theories  of  BaTi03  unify  Devonshire's  theory  and  Slater's  theory  by 
taking  into  account  the  contributions  by  all  ion  motions  to  spontaneous  polarization 

2  4-  I 

[Ref.  4],  As  we  see  from  Figure  2.4,  Ba.  and  Ti  ions  are  displaced  relative  to  the 
O  ions,  wrhere  produces  a  dipole  moment  and  then  spontaneous  polarization.  The 
upper  and  lower  oxygen  ions  may  move  downward  slightly  [Ref.  4]. 


Figure  2.3  Ionic  motions  in  BaTiO^  [Ref.3]. 

The  temperature  dependence  of  spontaneous  polarization  was  measured  by  Merz 
[Ref.  7],  from  hysteresis  loops.  This  is  illustrated  in  Figure  2.4. 

2.  Dielectric  constant  of  BaTiO^ 

As  is  seen  from  Figure  2.5,  the  tetragonal  phase  has  two  dielectric  constants 
that  are  denoted  by  ea  and  ec  .  This  is  a  result  of  the  crystal  symmetry  and  the  tensor 
character  of  the  dielectric  constant. 
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Figure  2.4  Temperature  dependence  of  spontaneous  polarization  [Ref.6]. 

In  the  cubic  phase,  there  is  only  one  dielectric  constant  and  above  the  Curie 
temperature,  the  dielectric  constant  obeys  the  Curie-Weiss  law  given  by  the  relation 

£  =  (  C/(T-Tc  )}.  (eqn  2.11) 

In  the  last  relation,  C  is  the  Curie  constant  ~  lCP  °K  and  T.  is  the  Curie  temperature. 
As  reported  by  Merz,  dielectric  constant  measurements  can  be  made  precisely  on  single 
domain  crystals.  At  low  temperatures,  the  orthorhombic  and  rhombohedral  phases  of 
BaTiO-  have  two  dielectric  constants  corresponding  to  the  symmetry  of  the  crystal. 

3.  General  ferroelectric  properties  of  BaTiO-j 

This  section  summarizes  the  ferroelectric  properties  of  BaTiO^  1)  It  has 
permanent  spontaneous  polarization  below  the  Curie  point  and  the  direction  changes 
with  each  phase.  The  magnitude  of  the  spontaneous  polarization  changes  with 
temperature  as  is  seen  from  Figure  2.6.  The  spontaneous  polarization  demonstrates 
hvsterises.  2)  As  we  see  from  Figure  2.6,  the  magnitude  of  the  spontaneous 
polarization  decreases  above  80°C.  At  the  Curie  point,  spontaneous  polarization 
vanishes.  3)  Spontaneous  polarization  is  in  the  direction  of  the  c-axis  of  the  crystal 
and  it  does  not  have  components  along  the  a-axis  of  the  crystal.  4)  The  dielectric 
constant  has  a  peak  at  the  Curie  temperature.  5)  Above  the  Curie  point  the  dielectric 
constant  falls  ofF  and  obeys  the  Curie-Weiss  law.  6)  At  the  Curie  point  a  ferroelectric 
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Figure  2.5  Dielectric  constant  of  BaTiO-j  as  a  function  of  temperature  T  [Ref.2]. 

to  paraelectric  phase  transition  occurs.  This  is  the  so-called  polarization  catastrophe 
[Ref.  4], 


Figure  2.6  Spontaneous  polarization  vs  temperature  [Ref.6]. 
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4.  The  pyroelectric  effect 

There  is  a  direct  relation  between  the  spontaneous  polarization  (key 
component  of  ferroelectricty)  and  the  pyroelectric  effect.  Lord  Kelvin  [Ref.  4]  defined 
the  pyroelectric  effect  as  follows:  "The  pyroelectric  effect  is  the  change  of  spontaneous 
polarization  due  to  temperature  changes".  So,  it  can  be  said  that  all  ferroelectrics  are 
pyroelectric,  but  the  reverse  is  not  true,  because  some  pyroelectric  crystals  may  not 
change  their  polarization  direction,  when  the  polarity  of  the  field  is  reversed.  The 
pyroelectric  effect  in  BaTiO-j  has  been  studied  by  Chynoweth  [Ref.  8].  When  the 
crystal  is  illuminated  ,  it  causes  a  small  change  in  temperature  of  the  crystal  which  in 
turn  causes  the  spontaneous  polarization  to  change.  In  other  words,  radiation  on  the 
crystal  causes  first  a  temperature  rise,  then  expansion  of  the  crystal  spacing,  and  finally 
a  change  in  spontaneous  polarization. 

When  electrodes  are  attached  to  the  polar  faces  (1  to  c-axis),  the  polarization 
produces  a  change  in  the  charge  on  the  electrodes,  explained  as  follows:  The  change  in 
the  charge  in  time  will  cause  the  current  to  flow  and  can  be  detected  by  an  external 
circuit.  To  develop  a  mathematical  relationship  consider  a  small  temperature  change 
dT  that  produces  a  change  dP$  in  P  .  If  this  temperature  change  occurs  in  a  time  dt, 
then  the  spontaneous  polarization  changes  at  the  rate  of  (dPs  /  dt),  which  is  equivalent 
to  a  current  density  J detected  by  the  external  circuit, 

Jpy  =  (dP$  /  dt)  =  (dPs/dT)  (dT/dt).  (eqn  2.12) 

For  a  small  temperature  change,  the  term  (dP$  /dT)  can  be  taken  as  a  constant  and  is 
called  the  "Pyroelectric  coefficient"  for  a  certain  temperature. 

p(T)  =  (dP$/dT)  (eqn  2.13) 

where  p(T)  is  the  pyroelectric  coefficient  [Ref.  9].  This  coefficient  can  be  found  for  a 
specific  temperature,  from  the  slope  of  the  P$  versus  T  curve,  at  that  temperature,  as  is 
seen  from  Figure  2.7. 

Since  p(T)  is  approximately  constant  for  a  specific  temperature,  the  current 
density  j7  can  be  detected,  if  (dT/dt)  is  large  enough.  We  know  that  the  spontaneous 
polarization  disappears  at  the  Curie  point  and  the  crystal  becomes  paraelectric.  Thus, 
the  pyroelectric  effect  should  not  be  observed  above  the  Curie  point,  but  has  been 
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observed  by  Chynoweth  [Ref.  10]  in  BaTiOj  .  He  concluded  that  BaTiO^  has  surface 
space-charge  layers  and  these  charges  produce  a  field  through  the  interior  of  the 
crystal.  Since  temperatures  above  the  Curie  point  were  not  used  in  this  research,  this  is 
beyond  the  scope  of  our  analysis.  The  derivation  of  the  Equation  2.12  will  be  given  in 
Section  H. 

C.  BULK  PHOTOVOLTAIC  EFFECT 

Chynoweth  [Ref  10]  has  reported  observing  the  photovoltaic  effect  in  BaTiO-j 
above  the  Curie  temperature.  It  has  been  interpreted  as  a  result  of  a  space-charge  field 
due  to  surface  charge  layers.  This  argument  led  Chen  to  explain  photorefraction  in 
electrooptic  crystals  incorrectly.  [Ref.  1].  Glass  [Ref.  11]  explained  the  microscopic 
origin  of  the  photovoltaic  effect  and  its  contribution  to  photorefraction.  Later,  Gower 
[Ref.  12]  observed  the  bulk  photovoltaic  effect  in  BaTiO-j  at  room  temperatures. 
Summarizing  the  previous  works,  the  observed  properties  of  bulk  photovoltaic  effect 
(BPVE)  in  BaTiO^  are  :  1)  Brody  [Ref.  13]  observed  that  the  photovoltage  (~ 

100-200V)  is  linearly  dependent  on  the  spontaneous  polarization  and  independent  of 
laser  intensity.  2)  Brody  [Ref.  13]  observed  that  the  photovoltage  vanishes  above  the 
Curie  temperature  (recall  spontaneous  polarization  vanishes  above  118°C  ).  3)  Koch  et 
al.  [Ref.  14]  have  observed  that  the  short  circuit  photovoltaic  current  (SCPVC) 
vanishes  at  temperatures  above  the  Curie  point,  and  the  SCPVC  is  linearly  dependent 
on  intensity.  4)  the  BPVE  exists  even  in  the  absence  of  an  external  electric  field  in  the 
pyroelectric  crystals.  Thermal  excitation  does  not  produce  the  BPVE. 

Now  with  the  help  of  these  observations,  we  can  compare  the  BPVE  with 
conventional  photovoltaic  effects.  In  general  there  are  two  types  of  photovoltaic 
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Figure  2.8  Photovoltage  vs  temperature  [Ref.  13]. 

effects.  The  first  one  is  known  as  the  Dember  or  photon  drag  effect  and  appears  when 
incident  radiation  on  the  surface  is  strongly  absorbed  by  the  semiconductor.  A 
potential  difference  developes  in  the  direction  of  the  radition  because  of  a  difference 
between  electron  and  hole  mobilities  [Ref.  15]. 

The  second  type  of  photovoltaic  effect  results  from  homogeneous  absorption  in 
the  crystal  containing  macroscopic  inhomogeneities,  as  in  semiconductor  p-n  junctions. 
Differences  in  the  Fermi  levels  bend  the  conduction  and  valence  bands. 

The  photovoltaic  effect  in  BaTiC>3  occurs  only  along  the  c-axis,  which  is  the 
spontaneous  polarization  direction.  In  our  experimental  geometry,  the  incident  beam 
direction  was  perpendicular  to  the  c-axis.  Thus,  the  possibility  of  the  Dember  effect 
was  ruled  out.  When  the  electrodes  were  attached  to  the  polar  faces,  which  were 
perpendicular  to  the  c-axis,  the  photo-voltage  across  the  polar  faces  was  around  10V. 
This  number  is  too  big  to  be  due  to  a  difference  in  Fermi  levels  of  a  p-n  junction 
mechanism  (around  IV).  A  possible  explanation  might  be  a  number  of  p-n  junctions 
connected  serially  to  give  10V.  The  macroscopic  homogeneity  of  BaTiCK  does  not 
allow  us  to  construct  a  model  of  this  kind;  so  this  possibilty  is  ruled  out. 

Thus,  it  was  seen  that  there  is  no  direct  anology  between  BPVE  and  the 
conventional  photovoltaic  effects.  So,  what  is  the  BPVE  and  what  kind  of  mechanism 
causes  this  effect  ?  In  order  to  answer  these  questions,  we  will  follow  Glass  et  al. 
[Ref.  16]  for  the  intrinsic  (pure)  crystal,  and  then  for  the  extrinsic  (impure)  case. 
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INTRINSIC  BPVE  :  Glass  et  al.  [Ref.  16]  considered  a  one  dimensional  polar 
structure,  in  which  a  unit  cell  contains  an  anion  O  and  a  cation  B.  In  BaTiO^  anions 
are  0“‘  ions  and  cations  are  Ti^+  ions.  The  top  of  the  valence  band  contains  oxygen 
2p  states  and  the  conduction  band  contains  Ti  3d  states.  BPVE  in  pure  (no  impurity) 
pyroelectrics  occurs  as  follows:  optical  transitions  are  delocalized  transitions,  provided 
that  there  is  an  electron-phonon  coupling  in  the  polar  crystal.  As  is  shown  in  Figure 
2.9,  GS  represents  the  ground  state  of  a  one  dimensional  polar  structure.  Consider  an 
optically  stimulated  excitation  that  causes  a  transfer  of  charge  from  an  anion  to  the 
nearest  cation.  Glass  et  al.  [Ref.  16]  assumed  that  a  charge  e  (<  <  electronic  charge) 
per  unit  cell  distributed  over  the  entire  chain  was  transferred  in  each  unit  cell.  In  the 
one  dimensional  argument,  charge  transfer  occurs  only  along  the  c  axis.  Figure  2.9b 
and  c,  shows  the  charge  distribution  after  excitation  in  the  Frank-Condon  state  (FCS) 
and  relaxed  state  (RES).  The  Franck-Condon  principle  says  that  when  an  charge  is 
excited  optically,  the  nuclei  of  the  ions  may  be  considered  to  remain  at  rest  during  this 
process.  In  other  words,  an  optical  excitation  occurs  in  a  very  short  time  compared 
with  the  period  of  corresponding  lattice  vibrations.  In  Figure  2.9c,  the  ionic  shift  has 
been  greatly  exaggerated. 
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Figure  2.9  shows  that  charge  moves  continuously  to  the  right  but  in  general,  the 
probability  of  charge  transfer  to  the  right  and  the  left  will  be  different  during  excitation 
and  recombination.  This  will  result  in  a  steady  state  photo  current.  This  model  looks 
like  a  microscopic  shift  register,  in  which  charge  £  is  transferred  from  one  unit  cell  to 
the  next  by  ionic  relaxation  in  a  random  manner.  Glass  et  al.  [Ref.  16]  have  performed 
calculations  based  on  these  ideas. 

EXTRINSIC  BPVE  :  This  case  is  similar  to  the  previous  one,  except  there  are 
impurities  in  the  crystal.  The  type  of  impurities  in  BaTiO^  will  be  mentioned  in  the 
next  section.  Baltz  et  al.  [Ref.  17]  took  up  this  case,  in  order  to  a  derive  short  circuit 
photovoltaic  current  (SCPVC)  quantum  mechanically,  which  was  originally  given  by 
Glass.  For  convenience,  we  will  follow  their  arguments  to  summarize  extrinsic  BPVE. 
They  have  considered  a  one  dimensional  model  of  a  photo-excited  carrier  in  a 
nonsymmetric  potential  V(x)  z  V(-x)  which  simulates  the  existence  of  an  acentric 
impurity.  In  Glass's  model,  the  SCPVC  is  due  to  electronic  transitions  from  a  localized 
ground  state  |0>  into  the  delocalized  states  |+  >  and  |->  ,  in  which  the  photo-excited 
carrier  moves  along  the  -l-c  and  -c  directions.  Since  the  potential  is  asymmetric,  the 
photoionization  probabilities  p+  and  p_  are  different.  This  results  in  currents  J+  in 
the  -He  direction  and  J_  in  the  -c  directions  subject  to 

(dJ/dt)  -H  (  1/t)J  =  (ea  I/hv  )vp  .  (eqn  2.14) 

This  equation  is  valid  for  J+  and  J  .  Here,  t  stands  for  the  life  time  of  the  carrier 
which  depends  on  the  recombination  and  collision  rates,  I  is  the  intensity  of  the 
incoming  beam,  a  is  the  absortion  coefficient  and  v  is  the  carrier  velocity.  In  the  steady 
state,  the  net  SCPVC  is  given  by 

Jph  =J+  -  J.  =(a  Ie/hv  )(1+  p+  -I.  p.  ).  (eqn  2.15) 

In  this  equation,  1+  =v+  t+  ,  1_  =v_  t_  and  1+  and  1_  are  the  mean  free  path  lengths 
of  the  excited  carriers.  In  the  case  of  thermal  excitation,  Jph  =0  ,  because  the 
probabilities  of  excitation  and  recombination  in  the  directions  -He  and  -c  of  the 
crysatal  are  the  same, 


p+  ( thermal)  =p_  (thermal)  =  exp(-E/kT). 


(eqn  2. 16) 


Here,  k  is  Boltzman's  constant  and  E  is  the  energy  of  the  particular  state.  In  the 
literature,  Jph  usually  is  given  by  this  form  [Ref.  1 1] 

Jph  =K(V  )a  (V)L  (eqn  2.17) 

In  the  last  equation,  K  is  the  anisotropy  constant  which  depends  on  the  frequency  v 
and  the  nature  of  the  absorption  centre.  It  represents  the  difference  in  probabilities 
times  the  mean  free  path  length  charge  moving  in  the  +  and  •  c  axis  direction.  It  is 
independent  of  the  crystal  geometry',  the  electrode  configuration,  and  the  impurity 
concentration.  Glass's  model,  predicts  the  mean  displacement  of  the  photo  carriers  per 
excitation,  which  is  restricted  to  the  lattice  constant. 


23 


D.  SEMICONDUCTOR  PROPERTIES  OF  BATIO3 

BaTiOj  has  a  very  small  conductivity.  Gower  [Ref.  12]  has  Found  a  conductivity 
of  ~  0.8*10*^  I  (ohm  cm)'*  ,  where  I  is  the  laser  intensity.  The  energy  band  gap  is  ~ 
3eV.  These  numbers  show  that  BaTiO^  is  a  semiconductor.  The  semiconducting 
properties  of  BaTiO^  are  thought  to  be  important  in  the  origin  of  photorefraction.  As 
has  been  reported  by  Klein  et  al.  [Ref.  18]  and  Feinberg  et  al.  [Ref.  19],  impurities 
determine  the  semiconductor  properties  of  BaTiO^  .  At  the  present  time,  however, 
little  is  known  about  the  impurities  that  are  responsible  for  photorefraction.  We  will 
summarize  the  experimental  results  of  Klein  and  Feinberg  about  the  impurity 
properties,  Before  geting  into  the  details,  it  will  be  useful  to  enumerate  some  basic 
definitions  for  semiconductors.  Donors  are  the  type  of  impurities  which  can  give  up 
electrons.  In  this  process,  since  the  electrons  still  stay  in  the  crystal,  the  entire  crystal 
stays  neutral.  Acceptors,  can  receive  electrons  from  the  valence  band,  in  order  to 
complete  the  covalent  bonds  between  neighboring  atoms.  In  this  case,  a  vacancy  or 


(  a  )7/////con ^LW//////////, 


A 

aonors  -0- 

-0-  0 

^  0- 
occepior* 

sue* 

+  ^  0 

note 

'///////;, 7; 

///.  voience  Dond  * / / / / , 

/ / / / / .  , .  t  / / / //// / 

r  L  - — ’ 

1 V////-/ 1 ' ‘ '  t  * '  f ///////////// 

'///  ,/ conduction  Don  a  ////  /////// 

a  on  or  1  -0 

-0"  /0 

1  >  4 

site*  / 

eteciron 

occepior* 

♦ 

7////// //////// ^77^7//// 


Figure  2.10  Donors- Acceptors  in  a  ferroelectric  semiconductor  [Ref.  19]. 
(a)  Majority  carriers  are  holes,  (b)  Majority  carriers  are  electrons. 
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hole  is  created  for  each  electron  removed  from  the  valence  band.  A  hole  acts  like  a 
positive  charge  under  an  applied  field.  The  mobility  is  the  magnitude  of  the  drift 
velocity  per  unit  electric  field  and  is  positive  for  electrons  and  holes,  even  though  their 
drift  velocity  is  opposite.  The  hole  mobility  is  smaller  than  the  electron  mobility 
[Ref.  4].  In  a  P  type  crystal,  acceptors  are  dominant  which  means  holes  are  the 
majority  carriers,  and  the  electrons  are  minority  carriers,  this  is  shown  in  Figure  2.10a. 
In  an  N  type  crystal,  donors  are  dominant,  so  the  electrons  are  majority  carriers  and 
the  holes  are  minority  carriers.  This  can  be  seen  from  Figure  2.10b. 

Klein  et  al.  [Ref.  18]  reported  three  major  groups  of  impurities  in  BaTiO^  :  1)  Ca 
and  Sr  are  in  the  same  group  with  Ba.  They  are  always  together  with  Ba  in  nature  and 
it  is  difficult  to  separate  them.  They  assumed  that  Ba  does  not  play  an  important  role 
in  photorefraction.  Therefore,  Ca  and  Sr  do  not  introduce  any  energy  levels  in  the 
energy  band  gap  of  BaTiO^  .  2)  Al  and  Si  exist  in  the  crystal  because  of  the  furnace 
walls  and  heating  elements  used  during  the  crystal  growth  process.  Since  they  have 
only  one  stable  valence  state  namely  Al3  +  and  Si^+  ,  they  can  not  be  involved  in  the 
intervalence  charge  transfer.  The  existence  of  transition  metals  in  many  oxide  crystal  is 
unavoidable  because  of  their  abundance  in  nature,  and  is  important  in  photorefraction. 
The  first  six  transition  metals  are  Cr,  Mn,  Fe,  Co,  Ni  and  Cu.  They  are  more  than 
likely  to  be  present  in  BaTiO^  .  Klein  et  al.  [Ref.  18],  have  done  absorption  coefficient 
measurements,  electron  paramagnetic  resonance  (EPR),  beam  coupling  and  spark 
source  emission  spectroscopy  (SIMS)  experiments.  These  experiments  showed  that  Fe 
is  the  impurity  responsible  for  photorefraction,  and  the  holes  are  the  majority  carriers. 
They  concluded  that  the  filled  acceptors  are  Fe3  +  and  the  empty  acceptors  are  Fe2  + 
ions.  The  ion  concentration  were  given  as 
N  =  (Fe2+  }  -  (2-9)*1016  cm'3  , 
and 

N+  =  (Fe3+  }  ~  (2-8)*  10 1 8  cm'3  . 

According  to  Klein  et  al.,  the  Fe  /  Fe  ratio  governs  the  photorefraction 
properties  of  BaTiO^  .  We  see  from  Figure  2.11  that  the  relative  contribution  of  holes 
and  electrons  to  photorefraction  depends  on  the  Fe  +  2  /Fe"1’3  ratio. 

An  alternative  interpretation  was  given  by  Feinberg  et  al.  [Ref.  19].  They  were 
able  to  alter  the  photorefractive  properties  of  a  nominally  pure  single  crystal  of  BaTiO-j 
by  treating  the  crystal  at  650°C  in  oxygen  at  different  partial  pressures.  This  treatment 
caused  the  effective  density  of  photorefractive  charge  carriers  to  change  in  the  crystal 
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Figure  2.1 1  Relative  contribution  of  electrons  and  holes  to  the  photoconductivity  [Ref.  18]. 

and  could  convert  a  dead  (inactive  in  photorefraction)  crystal  into  a  live  (active)  one. 
Reduction  (treatment  at  low  temperature)  decreases  the  Curie  temperature  of  the 
crystal  and  also  decreases  the  energy  band  gap,  implying  that  oxygen  vacancies  were 
added  to  the  bulk  crystal.  Negative  photorefractive  charge  donors  are  associated  with 
these  oxygen  vacancies,  depending  on  whether  the  partial  pressure  of  oxygen  was 
greater  than  or  less  than  0.5  atm  during  the  treatment.  Either  holes  or  electrons  may  be 
the  majority  carriers.  According  to  Feinberg  et  al.  [Ref.  19]  electron  and  hole 
contributions  to  the  photoconductivity  were  found  to  vary  with  the  density  of  oxygen 
defects  (vacancies)  introduced  into  the  crystal  during  the  treatment.  As  we  see  from 
Figure  2.12,  the  crystal  switches  from  N  type  to  P  type  at  an  oxygen  pressure  of  —  0.5 
atm  in  which  the  effective  density  N  of  the  photorefractive  charges  is  minimum. 
They  reported  that,  in  their  samples,  the  effective  density  of  photorefractive  charge 
carriers  was 

Npr  ~  (0.3-6.2)*1016  cm'3  . 

Between  these  two  studies,  there  is  a  major  discrepancy.  For  instance,  Klein  et  al. 
[Ref.  18]  assumed  that  the  oxygen  vacancy  can  trap  one  or  two  electrons,  leaving 
donor  levels  in  the  BaTiO-j  energy  band  gap.  But,  for  two  reasons,  these  levels  can  not 
participate  in  the  photorefraction:  1)  The  donor  levels,  which  are  occupied  by  oxygen 
vacancies  and  are  relatively  shallow,  can  be  thermally  ionized  at  room  temperature.  2) 
According  to  Klein  et  al.'s  beam-coupling  experiments,  the  holes  are  majority  carriers 
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TREATMENT  OXYGEN  PRESSURE  loim) 


Figure  2.12  Oxygen  vacancies  as  photorefractive  centres  [Ref.  1 9]. 

(a)  Density  of  photorefractive  carriers  vs  oxygen  pressure,  (b)  Relative  contribution 
(R)  of  holes  and  electrons  to  photorefraction. 

in  contradiction  to  the  donor  character  of  the  oxygen  vacancies.  However,  according  to 
Feinberg  et  al.  oxidation  and  reduction  treatments  only  change  the  density  of  oxygen 
vacancies  (donors)  and  not  barium  vacancies  (acceptors).  On  the  other  hand,  if  Fe  is 
an  active  photorefractive  impurity,  then  Fe~+  ions  are  donors  and  Fe^+  ions  are 
acceptors,  and  the  oxidation-reduction  treatments  would  have  changed  the  densities  of 
both  donor  and  acceptor  levels.  Feinberg  et  al.  [Ref.  19],  suggested  that  annealing 
experiments  should  be  performed  at  very  low  oxygen  pressure.  According  to  Feinberg 
et  al.  (oxygen  vacancy  picture),  this  treatment  'would  increase  the  density  of  donors 
(oxygen  vacancies)  without  changing  the  density  of  acceptors  (barium  vacancies)  and 
would  create  an  effective  photorefractive  charge  density.  However,  according  to  Klein 
et  al.  (Fe  picture),  this  treatment  would  cause  all  the  acceptors  to  be  converted  into 
donors,  resulting  in  no  effective  photorefractive  charge  density. 

ENERGY  LEVELS  OF  TETRAGONAL  BATIO3  : 

In  Figure  2.13,  CB  is  the  conduction  band  and  VB  is  the  valence  band,  and  Oh  and  C4 
v  stand  for  cubic  and  tetragonal  BaTiO^  .  The  top  of  the  valence  band  is  composed  of 
oxygen  2p  states  and  the  conduction  band  is  Ti  3d  states  for  the  cubic  and  tetragonal 
symmetries.  On  the  other  hand,  the  energy  levels  of  the  Ba  ion  do  not  play  an 
important  role  in  photorefraction  (according  to  Klein  et  al.,)  and  are  not  illustrated.  In 
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Figure  2.13  Energy  levels  of  Bn+  Og  +  clusters  in  BaTiO^  [Ref.18]. 

the  cubic  symmetry,  the  valence  band  edge  decomposes  into  states  labeled  tj  and  tj  u 
and  the  bottom  of  the  conduction  band  is  labeled  as  t2  g  ,  as  is  shown  in  Figure  2.13a. 
A  phase  transition  from  cubic  to  tetragonal  symmetry  results  in  split  and  shifted  energy 
levels  (see  Figure  2.13a).  Fe”"1’  (FeOg  )  and  Fe^+  (FeOg  )  have  a  complex 
energy  level  structure  in  BaTiO}  .  However,  Fe""1"  and  FeJ  complexes  do  have 
localized  states  in  the  energy  band  gap.  The  spin-up  (|  )  states  (generally  occupied)  are 
donors  and  spin-down  (j  )  states  (generally  empty)  are  acceptors,  as  is  seen  from 
Figure  2.13  b  and  c.  Thus,  with  this  picture  heteronuclear  intervalence  transitions 
involving  2p-»  Fe  3d  and  Fe  3d-»  Ti  3d  transitions  with  photon  energies  less  than  the 
energy  band  gap  are  possible. 

Hellwarth  et  al.  [Ref.  20]  showed  that  electrons  and  holes  can  exist 
simultaneously  under  the  illumination  of  the  crystal. 

E.  CHARGE  TRANSPORT  IN  BATIO3 

Charge  transport  may  be  caused  by  drift,  diffusion  and  the  bulk  photovoltaic 
effect  (BPVE)  mechanisms.  Currently,  there  are  two  models  in  which  these  mechanisms 
are  used  in  different  ways  to  explain  photorefraction. 

BAND  MODEL:  This  model  has  been  proposed  by  Kukhtarev  et  al.  [Ref.  2]  and 
Young  et  al.  [Ref.  21].  In  the  previous  section,  the  simultaneous  existence  of  electrons 
and  holes  has  been  mentioned.  So,  charge  transport  is  possible  by  both  electrons  and 
holes.  The  band  model  however  is  based  on  only  electron  transport  or  only  hole 
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transport.  It  does  not  handle  electron  and  hole  transport  simultaneously.  According  to 
the  band  model,  electron  transport  is  a  three  step  process:  1)  Photoionization  of  an 
electron  from  filled  state  (donor)  to  the  conduction  band.  2)  Migration  of  the 
photoexcited  electrons  to  dark  regions  in  the  crystal  can  be  caused  by  drift  and 
diffusion  through  the  conduction  band.  3)  Finally,  recombination  at  an  empty  state 
(ionized  donor).  Hole  transport  is  also  a  three  step  process,  in  this  case  X4-  is  an 
acceptor  or  filled  state.  The  three  steps  are:  1)  Photoionization  of  a  hole  from  an 
acceptor  X"1"  to  the  valence  band.  2)  Migration  of  photoexcited  hole  to  dark  regions 
by  diffusion  and  drift,  through  the  valence  band.  3)  Recombination  at  an  ionized 
acceptor  or  empty  trap  site  X.  Figure  2.14  summarizes  these  ideas. 


CONDUCTION  BAND 


y 


h 


y 


VALENCE  BAND 


Figure  2.14  Band  model  of  charge  transport  in  BaTiOj  [Ref.  18]. 

According  to  Klein  et  al.  X  corresponds  to  Fe^+  and  X+  corresponds  to  Fe^’4’  .  On 
the  other  hand,  according  to  Feinberg  et  al.  X  corresponds  to  oxygen  vacancies  and 
X"1”  corresponds  to  barium  vacancies.  If  electrons  diffuse  away  from  the  illuminated 
region,  the  resulting  space  charge  field  will  be  directed  to  the  dark  regions,  as  is  seen 
from  Figure  2.15a.  If  holes  diffuse  away  from  the  illuminated  region,  then  the  space 
charge  field  will  be  directed  to  the  illuminated  regions. 

As  we  see  from  Figure  2.15a  and  2.15b,  the  sign  of  the  space-charge  field  is  opposite 
for  the  two  charge  carriers. 

HOPPING  MODEL:  This  model  has  been  proposed  by  Feinberg  et  al.  [Ref.  22]. 
According  to  this  model,  carriers  move  by  hopping  from  a  filled  state  to  a  neighbor 
empty  state  (trap).  This  model  assumes  single  hops  (photo  assisted  tunneling)  between 
filled  and  empty  sites.  The  semiconductor  can  conduct  in  the  impurity  band  (the 
energy  levels  of  the  impurities  in  the  energy  band  gap)  by  electrons  hopping  from  a 
filled  donor  to  an  ionized  donor.  Impurity  band  conduction  sets  in  at  lower  donor 
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Figure  2.15  Redistribution  of  charges  under  illumination  [Ref.  19]. 

density  levels  if  there  are  acceptors  present  to  ionize  some  of  the  donors.  Under 
illumination,  electrons  from  donors  will  prefer  to  hop  to  an  ionized  (empty)  donor, 
rather  than  to  a  filled  donor  atom. 

These  two  models  predict  the  same  results  as  long  as  the  electron  (hole) 
recombination  time  is  short  compared  to  any  optical  excitation,  and  the  relative  charge 
modulation  is  small.  The  band  transport  model  has  been  used  widely  in  the  literature. 
This  model  will  be  followed  from  now  on. 

DRIFT:  Photoexcited  carriers,  electrons  in  the  conduction  band  or  holes  in  the 
valence  band,  can  move  freely  under  an  external  electric  field.  Carriers  migrate  out  of 
the  illuminated  region,  until  they  are  trapped  by  empty  traps  near  the  beam  periphery. 
An  external  field  produces  drift  of  photoexcited  carriers  through  the  band.  Positive 
and  negative  charge  separation  during  drift  gives  rise  to  a  space-charge  field  and  this  in 
turn  acts  as  feedback  to  redistribute  the  charges.  The  drift  current  density  Jdr  is  given 
by 


Jdr  =  e^e  nE  +  e**h  PE  ' 


(eqn  2.18) 
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In  the  last  equation,  ne  .  f*tl  are  electron  and  hole  mobilities;  n,  p  are  electron  and  hole 
densities.  E  =  Eext  +  Esc  ,  where  E  is  the  external  field  and  Esc  is  the  photoinduced 
space-charge  field  (vector  signs  will  be  dropped  from  now  on,  considering  the  c-axis 
direction).  Drift  charge  transport  is  caused  by  photoconductivity  cr  , 

Jdr  =  <*  E  ,  <T  =efie  n  +enh  p.  (eqn  2.19) 

The  charge  transport  length  in  this  process  is,  Le  =|i  t  E,  where  t  is  lifetime  of  the 
carrier. 

DIFFUSION:  Due  to  thermal  motion  of  the  excited  carriers,  charge  will  diffuse 
depending  on  the  gradient  of  the  carrier  concentration. 

Jjif  =  e^e  (dn/dz)  -  eDh  (dp/dz)  (eqn  2.20) 

In  this  equation,  De,h  =  ne,h  (kT/e)  are  the  Einstein  diffusion  coefficients  for  electrons 
and  holes.  Photoexcited  charges  travel  a  characteristic  distance  which  is  called,  the 
diffusion  length, 

Ld  =  (D^t)1'2.  (eqn  2.21) 

In  the  steady  state,  the  space-charge  field  will  saturate,  satisfying  the  relation 

Jdr  +  jdif  -  0  .  (eqn  2.22) 

This  can  be  written  explicitly  in  the  form 

en(z)n  (z)  E(z)  +  eD(dn/dz)  =0  .  (eqn  2.23) 

With  uniform  illumination  and  zero  external  field  (Eext  =0  ;  E=E$c  )  and  with  the 
charge  density  along  the  z  axis(c-axis)  given  by  n(z)=sinKz,  where  K  is  a  wave 
number,  the  space-charge  field  is 

E$c  =(kT/e)(K  cosKz/sinKz).  (eqn  2.24) 
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The  last  equation  alone  can  not  cause  photorefraction,  because  BPVE  was  not  taken 
into  consideration. 

BPVE:  As  explained  before,  this  effect  was  described  by  the  photocurrent  density 
Jph  =  K  a  I  .  Now,  the  general  charge  transport  can  be  written  in  terms  of  a  total 
current  density  J  as 


J(z,t)=<r  E(z,t)  + eD{dn(z,t)  dt}  +  k  a  I(z,t). 


(eqn  2.25) 


In  this  equation,  the  first  term  is  due  to  drift,  the  second  term  is  due  to  diffusion  and 
last  term  is  due  to  BPVE.  Depending  on  the  magnitude  of  these  terms,  one  of  them  can 
dominate  the  others.  In  this  equation,  n(z,t)=nd  +n,  (z,t),  where  nd  is  the  thermally 
excited  carrier  concentration  in  the  dark  and  nj  (z.t)  the  excited  carrier  concentration  in 
the  illumination.  First  Young  et  al.  [Ref.  21]  and  later  Kukhtarev  et  al.  [Ref.  2]  have 
discussed  the  principles  of  photorefraction  theory.  The  common  concept  of  these 
studies  was  based  on  two  crucial  points:  1)  A  dynamic  theory  has  to  be  included  in  the 
model  in  order  to  explain  the  oscillatory  behavior.  2)  The  effect  of  the  space-charge 
field  feedback  on  the  redistribution  of  electrons  has  to  be  considered  to  predict  the 
saturation  of  the  field.  This  can  be  done  by  including  the  dark  conductivity  in  the 
model.  According  to  Kukhtarev  et  al.  electrons  are  the  majority  carriers  and  they  will 
be  excited  from  trap  centers  (filled  states)  to  the  conduction  band  when  the  crystal  is 
illuminated.  Electrons  in  the  conduction  band  will  migrate  out  of  the  illuminated 
regions,  under  drift,  diffusion  or  BPVE,  until  they  are  traped  by  empty  states  around 
the  beam  periphery.  The  rate  of  electron  excitation  into  the  conduction  band  is 
{sl(z,t)+p  }Nj  ,  where  s  is  the  cross  section  for  photoionization,  p  the  thermal 
excitation  rate  of  electrons,  Nd  the  concentration  of  filled  states  (donors).  The 
recombination  (traping)  rate  is  yr  nNa  ,  where  yr  is  the  recombination  constant,  n  is 
the  concentration  of  electrons  in  the  conduction  band,  N  is  the  concentration  of 
empty  states  (trapping  centres  or  ionized  donors).  Photoionization  of  the  donors  will 
cause  charge  transport  and  recombination.  Then,  the  total  change  of  n  will  be 
determined  by  the  continuity  and  rate  equations.  The  continuity  equation  is  given  by 


dn(z,t)/dt=  dN'a  (z,t)/dt-(l/e)(dJ(z,t)/dz). 


(eqn  2.26) 


The  rate  equation  is  given  by 
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dNa  (z,t)/dt  =  (sl(z,t)+  P)(N(z,t)-Na  (z,t))-yr  n(z,t)Na  (z,t). 


(eqn  2.27) 


In  the  last  equation,  N=Na  ■+"  N,  is  the  total  concentration  of  dopants,  and  it  is 
assumed  that  only  a  small  part  of  the  donors  will  be  ionized.  The  excited  carrier 
lifetime  is  T  ={l/(y  Na)}  and  is  assumed  to  be  independent  of  n.  In  thermal 
equilibrium,  the  thermal  excitation  rate  is  equal  to  the  recombination  rate  P  N=nj  t  . 
If  these  relations  are  invoked  in  the  rate  equation,  then  it  takes  the  form, 


dNa  (z,t)/dt  =  g(z)-(n(z,t)-nd)/t 


(eqn  2.28) 


In  the  last  equation,  g(z)  is  the  generation  rate,  which  is  assumed  to  be  proportional  to 
the  light  intensity  I.  The  generation  rate  is  usually  given  by 


g(z)  =  («  Q)/(hv)I(z). 


(eqn  2.29) 


In  the  last  equation,  Q  is  the  quantum  efficiency,  h  is  Planck's  constant  and  a  is  the 
absorption  coefficient.  Inserting  these  relations  in  the  continuity  equation  gives 


dn(z,t)/dt=g(z)-((n(z,t)-nd  )/t}-(l/e){d  J(z,t)/d  z). 


(eqn  2.30) 


In  the  last  equation,  the  first  term  represents  the  generation  of  photocarriers,  per  unit 
volume  per  second  due  to  the  photoexcitation  process.  The  second  term  represents  the 
number  of  carriers  trapped  per  second  per  unit  volume  due  to  the  recombination 
process.  The  last  term  stands  for  carriers  leaving  a  unit  volume  per  second,  due  to  a 
current  density.  At  any  point  and  time,  the  accumulation  rate  of  the  space-charge 
density  is  given  by  the  continuity  equation  in  one  dimension, 


d  psc  (z,t )/d  t  =  -d  J(z,t)/d  z. 


(eqn  2.31) 


Combining  the  last  relation  with  Poisson's  equation  (V  ■Esc=  psc  /£  )  gives 


d  E$c  /d  z=-(l/£  )JQ  1  ( d  J(z,t )/d  z)dt  +  F'(t). 


(eqn  2.32) 
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Integration  of  the  last  equation,  with  respect  to  z,  gives  the  photo-induced  space- 
charge  field, 

Esc  (z,t)  =  -( 1/c  )f0  1  J(z,t)dt  +  F(t).  •  (eqn  2.33) 

In  the  last  equation,  F(t)  can  be  determined  from  the  boundary  conditions,  as  in  the 
following  case.  A  constant  voltage  source,  of  voltage  V,  is  applied  across  the  crystal,  so 
that  E  =  V/L,  where  L  is  the  difference  between  electrodes.  Integration  of  the  total 
field,  E(z,t)  =  V  L  +  E$c  ,  must  still  give  V, 

J0  1  E(z,t)dz  =  V,  (eqn  2.34) 

so  that 


Jo  ‘  Esc  (z>t)dz  =  0. 


(eqn  2.35) 


This  relation  is  true,  even  in  the  short  circuit  configuration. 

Recently,  Young  et  al.  [Ref.  23]  found  expressions  for  E$c  ,  by  solving  Equation 
2.31  and  2.33.  Their  model  does  not  include  the  contribution  of  holes  in  the  general 
transport  equation.  On  the  other  hand  Valley  [Ref.  24]  almost  found  a  complete 
expression  for  E$c  .  This  model  does  not  include  BPVE,  so  that  it  can  not  explain  the 
steady-state  condition.  A  complete  dynamic  theory  has  to  be  able  to  solve  the 
continuity,  rate,  general  transport,  Poisson's  and  wave  equations,  simultaneously.  Also 
the  general  transport  equation  has  to  include  diffusion,  drift  ,  BPVE  processes  and  the 
contribution  of  holes. 

F.  ELECTROOPTIC  EFFECT 

The  electrooptic  effect  is  the  change  of  the  index  of  refraction,  due  to  an  external 
field  or  space  charge  field.  We  know  that  BaTiCK  is  a  ferroelectric  crystal  with  a  phase 
transition  at  the  Curie  point.  Between  118°C  and  5°C,  it  has  tetrogonal  symmetry 
which  is  noncentrosymmetric  with  point  group  4mm.  In  this  phase  the  linear  (Pockel) 
electrooptic  effect  is  important.  Above  the  Curie  point  the  crystal  becomes 
centrosymmetric  with  point  group  m3m  and  is  cubic.  In  this  phase  the  Pockel  effect 
disappears.  So,  it  is  clear  that  symmetry  plays  an  important  role  in  the  electrooptic 
effect. 
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Next  consider  the  electromagnetic  wave  propagation  in  the  crystal.  The  media 
through  which  the  electromagnetic  wave  propagates,  determines  the  characteristic 
behavior  of  the  wave.  The  media  we  are  dealing  with  in  this  thesis  is  a  macroscopically 
homogeneous,  anisotropic  and  nonlinear  crystal. 

Electromagnetic  wave  propagation  in  the  crystal  is  completely  described  by  the 
index  ellipsoid  or  indicatrix  [Ref.  25].  In  general,  the  indicatrix  is  given  by 

[X2  /  nj2  }  +  {Y2  /  n22  }  +  {Z2  /  n32}  =  1  .  (eqn  2.36) 

In  the  last  equation,  X=Dx  /(£q  Ex  Dx  )l'/2  ,  Y=D  /( £q  E  D  )^2  , 

l  ■'?  "  y  j 

z=d,  /(ca  e7  d  y  1  ,  where  D  ,  D  and  D,  are  the  components  of  the  displacement 

z  v  z  z  Ay  z 

vector,  and  nj  ,  n2  and  n3  are  the  refraction  indices  along  the  principal  axes  of  the 
ellipsoid.  In  the  presence  of  the  electric  field  in  the  crystal,  the  indicatrix  equation  is 
modified 

an  X2  +  a22  Y2  +  a33  Z2  +  2a23  YZ  +  2a31  ZX  +  2a12  XY=  1 

(eqn  2.37) 

In  this  equation,  a;k  are  the  polarization  constants.  The  linear  change  in  the 
polarization  constants  due  to  a  field  is  given  by 

aik  -  U/  n;  nk  }5ik  =  I,  rikl  E,  .  (eqn  2.38) 

In  the  last  equation,  the  summation  includes  1=1  to  1=3.  The  6;k  is  the  Kroneker 
delta  function,  rik)  is  the  electrooptic  tensor  and  n;  is  the  diagonal  element  of  the  index 
of  refraction  tensor.  Following  Johnston  [Ref.  26],  the  index  of  refraction  change 
along  the  c-axis  is 

8  (n3  )  =  [n33  r33  8  PJ  /  {2(  £33  -1)},  (eqn  2.39) 

and  the  index  of  refraction  change  along  the  a-axis  is 

8  (nj  )  =  {nj3  r,3  8  P3  }  /  {2(  £„  -  1)}.  (eqn  2.40) 
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In  these  equations,  £33  is  the  dielectric  constant  in  the  direction  of  the  c  axis;  and  r^ 
=  rji3  and  £33  =  r333  are  the  electrooptic  coefficients  along  the  a  axis  and  along  the  c 
axis,  respectively.  Johnston's  [Ref.  26]  model  does  not  fully  explain  photorefraction, 
because  of  the  following  problems:  1)  It  can  not  explain  the  existence  of  the  bulk 
photovoltaic  steady-state  current.  2)  According  to  this  model,  the  open  circuit 
photovoltage  should  decay  with  the  dielectric  relaxation  time  of  the  crystal.  In  reality, 
this  voltage  eventually  reaches  a  steady-state  value  (see  Figure  2.8).  3)  This  model 
requires  an  extremely  high  density  of  free  electrons  to  cause  photorefraction.  Overall, 
this  model  should  be  modified  by  taking  into  account  the  major  cause  of 
photorefraction,  which  we  already  know  is  the  photo-induced  space-charge  field  that 
resultes  from  the  charge  transport  mechanisms.  Following  Kukhtarev  et  al.  [Ref.  2], 
the  photo-induced  refractive  index  modulation  is 


S  (n3  )  =  (1.2)n33  r33  Esc  , 


(eqn  2.41) 


and 


5(nj  )  =  (l/2)nj3  r13  Esc. 


(eqn  2.42) 
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G.  PHOTO  REFRACTIVE  EFFECT  IN  BATIO3 

Several  mechanisms  that  result  in  photorefraction,  were  mentioned  in  the 
previous  sections.  In  this  section,  all  these  ideas  will  be  put  together,  in  order  to 
explain  photorefraction.  Since  beam-coupling  is  a  significant  manifestation  of 
photorefraction,  it  will  be  taken  as  an  example. 


Figure  2.16  Formation  of  interference  pattern  [Ref.27]. 

As  is  seen  from  Figure  2.16,  the  interaction  of  two  coherent  beams  (S:  signal  beam,  R: 
reference  beam)  crossing  at  20q  (refraction  angle)  creates  an  intensity  interference 
pattern  along  the  c-axis.  As  we  see  from  Figure  2.16  the  fringe  spacing  can  be  varied 
by  changing  the  incident  angle  0  .  The  wave  number  of  this  pattern  is  Kg  =  2ksin0  , 
where  k  is  the  wave  number  of  the  incident  beam  k=  2tt  /X  .  The  intensity  interference 
pattern  can  be  expressed  by 

I(z)  =  lQ(l  +  mcosK  z).  (eqn  2.43) 

In  this  equation,  Iq  is  the  total  intensity  of  each  of  the  two  waves, and  m  is  the 
modulation  ratio.  Under  the  absorption  of  light,  electrons  will  be  excited  from  filled 
donors  to  the  conduction  band  and  holes  will  be  excited  from  acceptors  to  the  valence 
band.  Photocarriers  will  have  the  same  perodicity  as  the  interference  pattern.  These 
photocarriers  will  migrate  to  the  dark  regions,  under  the  influence  of  diffusion,  drift 
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and  the  BPVE  mechanisms.  Then,  charges  recombine  with  empty  traps.  Charge 
transport  results  in  a  space-charge  grating.  This  periodic  space-charge  distribution  has 
the  same  spatial  frequency  as  the  intensity  interference  pattern.  But,  in  general,  the 
space-charge  grating  is  out  of  phase  with  the  interference  pattern. 


Figure  2.17  Photorefraction  due  to  beam-coupling  [Ref.  18]. 

This  photo-induced  space-charge  distribution  p  (z)  results  in  a  periodic  space-charge 
field  Esc  as  a  result  of  Poisson's  equation.  Young  et  al.  [Ref.  28]  found  an  expression 
for  the  photo-induced  space-charge  field  Esc  ,  without  including  the  contribution  of 
holes  to  the  charge  transport  processes.  For  the  case  in  which  diffusion  is  the  only 
transport  mechanism,  Esc  is  given  by 

Esc  (z,t)=  {(eg(z,t)tm)/e}{(KLd2  )/(l+(KLd)2  )1/2  }sinKz  .  (eqn  2.44) 
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For  the  case  of  drift  alone,  E$c  is  given  by 

E$c  (z,t)  =  -((eg(z,t)tm)/£}{Le  /(I  +(KLe  )2  )1/2  }cos(Kz  +  <pe  ).  (eqn  2.45) 

For  the  case  of  BPVE  alone,  E$c  can  be  expressed  by 

Esc  (z,t)=  ((eg(z,t)tm)/£}{Lph  /(l  +  (KLph  )2  )l>2  }cos(Kz-<pph  ).  (eqn  2.46) 

In  the  last  three  equations,  K  =  K  is  the  interference  pattern  wave  number,  g(z,t)  is  the 

© 

charge  generation  rate,  and  z  =£33  is  the  extraordinary  dielectric  constant.  Ld  ,  Le  and 
Lph  are  transport  lengths  in  diffusion,  drift  and  in  the  BPVE  processes,  respectively. 
The  tpe  and  cpph  are  phase  shifts  from  the  interference  pattern  due  to  drift  and  the 
BPVE  processes,  respectively.  The  relations  between  transport  lengths  and  phase  shift 
angles  can  be  written  in  this  fashion 

<Pph  =  tarfl  KLph  ’  =tan*1  KLe  •  (ecln  2-47) 

In  the  last  equation,  the  charge  transport  length  in  the  BPVE  prosess  is  given  by  Lph 
=  1+  p+  -  1.  p.  .  From  these  relations,  it  can  be  seen  that  phase  of  the  Esc  is  strongly 
dependent  on  the  transport  length  in  the  process.  If  diffusion  dominates  over  the  other 
mechanisms  in  the  general  charge  transport,  then  it  leads  to  a  phase  shift  of  90° 
between  the  space-charge  field  and  the  intensity  pattern.  This  can  be  seen  comparing 
Equation  2.43  with  Equation  2.44. 

Finally,  the  space-charge  field  modulates  the  index  of  refraction  via  the 
electrooptic  (Pockel)  effect.  But,  we  know  that  there  are  small  (negligible)  contributions 
of  polarization  changes  to  photorefraction,  given  by  Equation  2.39 

6  (n3  )  =  (l/2)n33  r33  {E$c  (z,t)  +  8  P3  /( E33  -1)}.  (eqn  2.48) 

Since  Egc  is  out  of  phase  with  I(z),the  refractive  index  change  will  also  be  out  of  phase. 
This  phase  shift  between  the  interference  pattern  and  refractive  index  pattern  (its 
magnitude  depends  on  the  charge  transport  mechanism),  causes  energy  tranfer  between 
the  two  beams.  This  is  the  so  called  beam-coupling.  All  these  ideas  can  be  summarized 
in  Figure  2.17  If  the  incident  beam  intensities  are  not  equal  during  the  beam-coupling, 
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energy  will  be  transferred  from  the  strong  beam  to  the  weak  beam.  Valley  [Ref.  24] 
assigns  the  direction  of  the  energy  transfer  according  to  the  type  of  carriers  involved  in 
the  charge  transport.  For  the  geometry  depicted  in  Figure  2.16,  if  the  energy  transfer  is 
from  the  reference  beam  to  the  signal  beam,  then  the  majority  carriers  are  holes.  But  if 
the  energy  transfer  is  reversed,  then  the  majority  carriers  are  electrons.  Beam-coupling 
and  energy  transfer  are  strongly  dependent  on  the  transport  length  of  the  charge 
transport  process. 

COUPLED  WAVE  THEORY:  This  is  another  way  to  explain  beam  coupling.  If 
the  Bragg  condition  is  satisfied,  the  hologram  grating  will  diffract  the  incident  beams. 
The  fraction  of  the  incident  light  which  is  diffracted  is  related  to  the  amplitude  of  the 
sinusoidally  varying  index  change  6  n.  Kogelnik  [Ref.  29]  has  given  a  relation  for  the 
diffraction  efficiency  r| 

r|  =  sin2  (Ml),  (eqn  2.49) 

where  M  =  7t  6  n  /X  cos0()  .  Here,  6  n  is  the  magnitude  of  the  index  of  refraction 
change,  1  is  the  interaction  length  of  the  two  beams,  X  is  the  wavelength  of  the  incident 
light  and  0q  is  the  angle  of  refraction.  The  interaction  length  is  comparable  to  the 
crystal  thickness  and  can  be  considered  as  the  thickness  of  the  hologram  grating. 
Amodei  et  al.,  [Ref.  30]  have  applied  coupled-wave  theory  to  the  beam-coupling 
problem.  Important  aspects  of  this  work  are  given  by 

(1^  /I$  )=cos2  (Ml)  +  A2  sin2  (Ml)-Asin(2Ml)sin(p  ,  (eqn  2.50) 

and 

(Is1  /Is  }  =  sin2  (Ml)  +  A2  cos2  (Ml)  + Asin(2Ml)sin(p  .  (eqn  2.51) 

In  the  last  two  equations,  I  ^  is  the  intensity  of  the  strong  beam  exiting  the  crystal, 
when  the  weak  beam  is  turned  on.  I$*  is  the  intensity  of  the  weak  beam  exiting  the 
crystal  when  the  strong  beam  is  turned  on.  Is  is  the  intensity  of  weak  beam  exiting  the 
crystal,  when  the  strong  beam  is  turned  off.  A  is  the  ratio  of  I$  to  If  ,  where  If  is  the 
intensity  of  strong  beam  exiting  the  crystal  when  the  weak  beam  is  turned  off.  If  we 
subtract  the  last  two  equations  from  each  other,  having  used  the  fact  that  1^  =  If  -A  I 
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and  1^  =  I$  +  A  I,  then  an  amount  of  energy  A  I  is  transferred  from  the  strong  beam 
to  the  weak  beam,  as  given  by  the  relation  [Ref.  30]: 

(1/A)-2F-1  =  (1-A“  )cos(2Ml)-2Asin(2Ml)sin<p  .  (eqn  2.52) 

In  this  equation,  F=A  I/Ig  is  called  the  fractional  energy  transfer,  and  the  relation 
between  the  fractional  energy  transfer  and  diffraction  efficiency  is  given  by 

(A  I/Is  }  =  2r|1/2  •  (eqn  2.53) 

H.  DERIVATION  OF  THE  QUANTUM  EFFICIENCY  OF  BATIO3 

In  this  section,  the  pyroelectric  effect,  bulk  photovoltaic  effect  and  results  of 
coupled  wave  theory  will  be  used  to  determine  the  quantum  efficiency  of  BaTiO}  . 
There  have  been  some  attempts  to  find  the  quantum  efficiency  of  ferroelectric 
semiconductors  in  the  literature.  Glass  et  al.  [Ref.  31]  found  a  quantum  efficiency  0.17 
for  LiNbCU  .  A  diffraction  efficiency  method  was  used  to  get  that  number.  As  they 
noted,  a  quantum  efficiency  of  0.17  gave  a  trapped  electron  density  that  was  three 
times  the  Fe  concentration  of  the  crystal  (impossible).  Later,  Young  et  al.  [Ref.  32] 
found  a  quantum  efficiency  of  ~  10'^  for  LiNbO^  .  They  first  estimated  the 
photovoltaic  transport  length,  and  using  Equation  2.15,  calculated  the  quantum 
efficiency. 

First,  using  electrostatics,  we  shall  derive  the  pyroelectric  current  equation. 


Figure  2.18  Polarization  surface  charges. 
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Figure  2.18  is  appropriate  for  a  single  domain  crystal.  This  means  the  volume 
polarization  charge  density  is  zero  (p  =  -V  .  P  =  0).  On  the  other  hand,  the 
polarization  charge  densities  on  the  surfaces  are  not  zero  but  are  given  by 

Ps.Nl  =  Pplandfs.ff2-Pp2.  (eqn  2.54) 

In  the  last  equation,  N j  and  N ->  are  unit  vectors  normal  to  the  surfaces;  while  pp| 
=  qpl  /  Asj  and  pp->  =  qp2  /A$2  are  P°lar'zation  surface  charge  densities.  In  our 
crystal  the  surface  areas  are  equal,  A$j  =  A  ->  =A$  From  the  relations  above,  the 
surface  polarization  charges  are 

qp  i  =  ps  As and  v  =  -ps  As  •  (e9n  2-55) 

When  the  crystal  was  illuminated,  the  crystal  warmed  up  and  as  a  result  expanded  the 
unit  cell.  This  expansion  caused  a  change  in  the  spontaneous  polarization.  If  the  rate 
of  change  in  temperature  was  large  enough,  a  pyroelectric  current, 

ipy  =  (dqp  /d t)  =  As  (dPs  /dt)  (eqn  2.56) 

occured.  The  change  in  P$  was  due  to  heating, 

(dP$  ,'dt)-(dPs  /dT)(dT/dt).  (eqn  2.57) 

Invoking  this  equation  in  Equation  2.56  gives 

ipy  =A$(dPs/dT)(dT/dt).  (eqn  2.58) 

Using  the  definition  in  Equation  2.13,  we  can  rewrite  the  last  equation  in  the  form 

ipy  =A$p(T)(dT/dt),  (eqn  2.59) 

where  p(T)  and  (dT/dt)  are  the  pyroelectric  coefficient  and  the  rate  of  temperature 
change,  respectively.  We  need  to  know  how  much  power  was  absorbed  to  produce  this 
current.  We  can  find  this  with  a  thermodynamic  relation.  Thermodynamically  the 
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absorbed  power  equals  the  heat  to  warm  up  the  crystal  minus  the  radiative  heat  loss. 
The  radiative  heat  loss  obeys  the  Wien  radiation  law  and  is  negligible  because  of  the 
low  power  range  in  the  experiments:  the  crystal  temperature  did  not  change 
significantly  with  respect  to  room  temperature  ( ~  9°C).  The  absorbed  power  per  cm 
(®)  can  be  written  as 

<t>  =JCp  p  Ab(dT/dt),  (eqn  2.60) 

where  Cp  is  the  specific  heat  of  the  crystal  at  constant  pressure,  Ab  is  the  beam  area,  J 
is  the  mechanical  equivalent  of  heat  (i.e.,  Joules  per  Calorie),  and  p  is  the  density  of 
BaTiOj  .  We  can  rewrite  the  last  equation,  combining  it  with  the  pyroelectric 
equation,  as 

=  ((Jcp  P)/p)(Ab  /As  )ipy  .  (eqn  2.61) 

The  last  equation  relates  the  absorbed  power  per  cm  (<I>  )  to  the  pyroelectric  current. 
The  next  step  was  to  find  the  number  of  absorbed  photons  per  second  per  cm  (N  ) 

N  =(<!>  /hv  )  sec*^  cm'*  .  (eqn  2.62) 

In  the  last  equation,  h  is  Planck's  constant,  and  v  is  the  frequency  of  the  incident 
beam.  In  order  to  find  the  quantum  efficiency,  we  need  to  know  the  number  of  photo- 
excited  carriers,  corresponding  to  N  .  During  our  photo-induced  current 
measurements,  an  external  field  was  not  applied  (see  Section  III-B).  This  means,  in  the 
general  charge  transport  Equation  2.25,  that  the  first  term  is  zero.  The  remaining  two 
terms  in  this  equation  are  due  to  diffusion  and  BPVE.  The  absorption  coefficient 
measurements  showed  that  the  impurities  are  anisotropically  distributed  in  the  crystal 
(see  section  III-A).  As  already  discussed  in  Section  II-C,  the  BPVE  is  due  to  the 
anisotropically  distributed  impurities  in  the  crystal.  In  fact  a  short  circuit  photovoltaic 
current  was  observed  in  the  steady-state  in  our  photocurrent  measurements  (see 
Section  III-B).  Moreover,  Kratzig  et  al.  [Ref.  33]  were  able  to  reduce  the 
photorefraction  (optical  damage)  by  compensating  for  the  BPVE.  This  indicates  that 
the  BPVE  will  be  the  key  component  in  the  quantum  efficiency.  If  we  rewrite  Equation 
2.15,  and  using  2.17,  then  the  photo-excited  carriers  per  cm  per  second  N.  is 
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Xc-(iph/LPh  <0<A„/As). 


(eqn  2.63) 


In  this  equation,  iph  is  the  SCPVC,  Lph  is  the  transport  length  in  the  photovoltaic 
process,  defined  under  equation  2.47,  and  e  is  the  electronic  charge.  The  quantum 
efficiency  is  given  by 


Q-<VN  P  >• 


(eqn  2.64) 


Rearranging  Equation  2.65  in  terms  of  Lph  we  have 

QLPh  =(‘ph  ^  eXhv  XAb  'As  ) cm-  (ecin  2-65) 

The  surface  area  of  our  crysatal  is  As  =0.25cm*'  .  The  diameter  of  the  beam,  defined 
by  the  e'2  value  of  its  intensity  is  1.3mm  (as  given  by  the  manufacturer).  This  gives  Afe 
=  0.0132cmi'  .  The  rest  of  the  parameters  involved  in  Equation  2.62  can  be  read  from 
Table  1.  Using  J  =  4.26  (Joule'calorie),  C  =0.12  (calorie/gram  °C),  p  =  6  (gram/cm2  ) 

O  J  ^ 

and  p=  2*10'°  (Coulomb  cm  ),  Equation  2.62  becomes 

cp  =0.0792*  10+ 8  i  ,  =a  W  Watt/cm.  (eqn  2.66) 

py  p  n 

In  the  last  equation  W  is  the  incident  beam  power  and  ap  is  the  photorefractive 
absorption  coefficient.  Rewriting  Equation  2.18  as 


iph  =  K  ap  w(As /Ab  ),  (eqn  2.67) 

and  using  Equation  2.67  we  find  the  anisotropy  constant  K  to  be 

K  =  (iph  /<P  )(Ab  / A$  )((Ampere-cm)/watt).  (eqn  2.68) 

If  we  rewrite  Equation  2.66,  it  becomes 

QLph  =  K  (v)  ((hv  )/e)  cm.  (eqn  2.69) 


44 


This  relation  was  used  by  Young  et  al.  [Ref.  34]  and  Vahey  [Ref.  35],  but  they  all  took 
large  values  for  the  absorption  coefficient  in  the  calculation  of  the  quantum  efficiency 
of  LiNbOj  .  To  the  best  of  our  knowledge,  this  is  the  first  reported  calculation  of  Q  for 
BaTiC>3  . 

The  next  task  was  to  find  the  carrier  transport  length  Lph  in  the  photovoltaic 
process.  The  results  of  coupled  wave  theory  were  used.  As  already  stated,  the  amount 
of  energy  transferred  from  the  strong  beam  to  the  weak  beam  depends  on  the  phase 
angle  between  the  hologram  grating  and  the  interference  pattern.  This  phase  angle  in 
turn  is  determined  by  the  carrier  transport  length  Lph  .  If  Equation  2.53  is  rearranged, 
the  phase  shift  angle  (p  is 

<p  =  sin'1  (( 1-A2  )cos(2Ml)-((  l/A)-2F-l)}/{2Asin(2Ml)}.  (eqn  2.70) 

Since  BPVE  dominates  in  the  general  charge  transport  equation  (due  to  Equation 
2.46),  (p  =<Pph  •  Arranging  Equation  2.47,  the  charge  transport  length  is 

LPh  a,0/K)(tanq>ph  )•  (efin  2-71) 

Combining  the  last  equation  with  Equation  2.70,  the  quantum  efficiency  is 

Q(v  )  =  (K  (v  )hv  K)/(etan(pph).  (eqn  2.72) 

In  the  last  equation,  Q  is  the  quantum  efficiency,  K  is  the  wave  number  of  the  grating, 
h  is  Planck's  constant,  v  is  the  frequency,  e  is  the  electronic  charge  and  (pph  is  the  angle 
betw’een  the  hologram  grating  and  the  intensity  pattern. 

In  order  to  find  an  actual  number  for  the  quantum  efficiency  Q,  we  need  to  have 
the  values  of  the  photorefractive  absorption  coefficient  ap  ,  the  short  circuit 
photovoltaic  current  iph  ,  the  pyroelectric  current  ipy  and  the  phase  shift  angle  between 
the  hologram  grating  and  the  intensity  grating  (pph  .  These  parameters  had  to  be 
obtained  from  experiments.  The  photocurrent  measurements  w’ere  made  to  determine 
ap  ,  ipy  and  iph  .  Beam-coupling  measurements  wrere  made  to  determine  (pph  .  In  the 
next  chapter,  the  experimental  methods  related  to  the  microscopic  parameters 
mentioned  above  will  be  presented.  At  the  end  of  the  next  chapter,  the  quantum 
efficiency  will  be  calculated,  using  the  theory  presented  here,  and  the  experimentally 
determined  parameters. 
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TABLE  1 

MICROSCOPIC  PARAMETERS  OF  BATIO3 


PARAMETER 

SYMBOL 

VALUE 

REFERENCES 

Size  of  tetragonal 

c 

4.036A0 

[Ref.2] 

unit  cell. 

a 

3.992A0 

" 

Dielectric 

ec 

160 

[Ref.  6] 

constant 

Ca 

4100 

Mag. of  spon. 
pol. vector. 

Ps  16 

*10‘6  Coulomb  cm‘“ 

" 

Mass  density 

P 

6  gr  cm*-* 

[Ref.9] 

Specific  Heat 

S 

0.12  cal/gr°C 

- 

Mec.Eq.of  Heat 

J 

4.16  Joule/calorie 

- 

Pyroelectric  coef. 

P  2 

*10*^  Coulomb  cm'" 

» 

Index  of 

no 

2.437 

[Ref.2] 

Refraction 

ne 

2.365 

- 

Electrooptic 

rll 

8  pm/V 

[Ref.35] 

r33 

23  pm/V 

coefficients 

r42 

820  pm/V 

Recomb. cons. 

Yr 

5*10*^  cm*  /sec 

[Ref.35] 

Mobility 

P 

0.5  cnV*  /V sec 

[Ref.  34] 
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III.  EXPERIMENTAL  OBSERVATIONS 


A.  ABSORTION  COEFFICIENT  MEASUREMENTS 

1 .  Purpose  of  this  experiment 

Absorption  coefficient  measurements  play  an  essential  role  in  identifying  the 
types  of  impurities  in  a  crystal.  However,  it  is  not  the  only  method  since  EPR 
measurements  can  achieve  for  the  same  purpose  (Ref.  18],  Our  interest  is 
photorefraction  and  the  subsequent  beam-coupling  at  visible  wavelengths.  Absorption 
coefficient  measurements  are  helpful  in  identifying  the  the  band  gap  and  the  type  of 
impurities  in  semiconductor  crystals.  If  one  had  access  to  a  crystal  doping  facility,  one 
could  dope  the  crystal  with  a  suspected  impurity  and  observe  the  effect  on  the 
absorption  spectra.  Unfortunately,  we  did  not  have  the  doping  facilities  to  do  this. 
Absorption  coefficient  measurements  helped  us  find  the  energy  band  gap  and  to 
understand  the  molecular  vibration  modes  in  BaTiO^  .  The  infrared  absorption 
spectrum  gave  information  on  the  vibrational  motion  of  the  Tf*+  -0“'  ions,  which  are 
the  ions  responsible  for  spontaneous  polarization.  Not  all  of  the  absorbed  energy 
contributes  to  photorefraction.  We  suspect  that  some  of  the  absorbed  energy  is  lost  on 
the  internal  reflections,  scattering  and  etc. 

2.  Procedure 

Ultraviolet  (UV),  visible  and  near  infrared  (NTR)  absorption  coefficient 
measurements  were  made  by  using  a  Perkin-Elmer  Model  330  spectrophotometer.  It 
covers  the  187-2500nm  wavelength  range  and  the  wavelength  accuracy  is  0.2nm  for  the 
UV  region  and  Inm  for  the  NIR  region.  The  output  light  beam  of  this  instrument  was 
partially  polarized.  For  our  purposes,  polarized  light  was  needed.  Therefore,  a  polarizer 
was  used  to  get  horizontally  polarized  light.  The  spectrophotometer  beam  size  was 
larger  than  the  crystal  size,  so  a  special  holder  was  used  to  confine  the  light  beam  to 
the  crystal  size. 

In  order  to  measure  the  absorption  coefficient  of  the  crystal,  the  absorption 
due  to  the  polarizer  and  holder  had  to  be  eliminated.  This  was  achieved  by  simply 
making  a  correction  run  before  the  measurement.  Therefore,  without  the  crystal,  the 
instrument  indicated  zero  absorption. 
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In  the  infrared  (IR)  region  of  the  spectrum,  absorption  coefficient 
measurements  were  made  by  using  a  Perkin-Elmer  Model  337  IR  spectrophotometer. 
This  has  a  2.5-25p  m  wavelength  range.  In  this  case,  a  polarizer  was  not  used  but  a 
beam  restricting  aperture  was  used.  This  instrument  was  not  computerized,  so  a 
correction  run  was  not  needed.  Instead,  a  shutter  was  used  on  the  reference  beam  path, 
and  then  the  gain  of  the  instrument  was  adjusted  to  give  zero  absorption  without  the 
crystal  in  the  signal  beam  path.  These  spectrophotometers  measure  the  absorbance, 
which  is  dimensonless.  In  order  to  get  a  absorption  coefficient  a  ,  this  relation  has  to 
be  used  «  =  Absorbance  crystal  length  .  Our  crystal  is  0.5  cm  long.  Then  the 
absorption  coefficient  a  =  (2*Absorbance)cm'^  . 

The  reflection  loss  from  a  surface  of  index  of  refraction  n  is  given  by 

R  =  (n-1)2 /(n+1)2  .  (eqn  3.1) 

In  this  equation  R  is  the  ratio  of  the  reflected  beam  intensity  to  the  incident  beam 
intensity  and  n  is  the  index  of  refraction  at  a  wavelength  where  the  crystal  is 
nonabsorbing.  R  does  not  vary  appreciably  with  wavelength.  Using  n  =  2.4,  R  was 
found  to  be  0.17  [Ref.  36].  In  the  calculation  of  quantum  efficiency,  the  absortion 
coefficient  obtained  from  this  method  were  not  used,  instead  we  used  absorption 
coefficients  derived  the  pyroelectric  measurements.  Because  of  this,  Fresnel  reflection 
losses  were  not  removed  from  data.  These  measurements  were  made  at  a  room 
temperature  of  23°C. 
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3.  Results  of  the  experiments 


Figure  3.1  Measured  absorption  spectrum  of  BaTiO^  ,  E  . L  c,  K  //  c. 
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Figure  3.2  Measured  absorption  spectrum  of  BaTiO-j  ,  E  1  c,  K  //  c. 
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Figure  3.3  Measured  absorption  spectrum  of  BaTiO^  ,  E  1  c,  K  1  c. 
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Figure  3.4  Measured  absorption  spectrum  of  BaTiO^  ,  E  1  c,  K  J.  c. 
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Figure  3.5  Measured  absorption  spectrum  of  BaTiOj  ,  E  //  c,  K  1  c. 
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Figure  3.6  Measured  absorption  spectrum  of  BaTiO^  ,  E  //  c,  K  1  c. 
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Figure  3.7  Measured  absorption  spectrum  of  BaTiO-j  at  IR  region. 
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4.  Interpretation  of  the  results 

As  we  see  from  Figure  2.5,  the  dielectric  constant  of  BaTiCA  is  anisotropic 
in  other  words,  £a  is  not  equal  to  £c  .  Because  of  the  tetragonal  crystal  symmetry  we 
can  expect  anisotropic  absorption  coefficients  for  BaTiO-j  .  In  the  370-900nm 
wavelength  range,  the  absorption  coefficient  measurements  were  made  at  six  different 
crystal  orientations.  The  electric  field  of  the  incident  beam  can  be  either  parallel  or 
perpendicular  to  the  c  axis  of  the  crystal  in  these  configurations.  For  instance,  in 
Figure  3.1  and  3.2,  E  1  c  and  K  //  c  ,  where  K  is  the  propagation  direction,  but  the  c 
axis  direction  was  reversed  in  Figure  3.2.  In  Figure  3.3  and  3.4,  E  1  c,  K  1  c,  but  the  c 
axis  direction  was  reversed  in  the  second  one.  In  Figure  3.5  and  3.6,  E  II  c  and  |  1  c 
but  in  Figure  3.6,  the  c  axis  direction  was  reversed.  For  each  configuration,  the 
absorption  spectrum  was  different.  This  is  an  indication  of  the  anisotropic  absorption 
center  (impurity)  distribution. 

Conurion  characteristics  of  Figure  3. 1-3.6  are:  1)  There  is  a  broad  absortion 

plateau  between  880nm  and  the  band  edge  at  460nm.  2)  The  maximum  absorption  is 

at  ~  410nm,  which  gives  us  the  band  gap  energy.  Feinberg  et  al.  [Ref.  19]  showed 

that  the  absorption  plateau  became  smaller  after  the  crystal  had  been  treated  with  0.2 

atm  of  0->  .  The  band-edge  shifted  to  shorter  wavelengths,  as  the  oxygen  pressure  of 

the  process  was  increased,  indicating  that  the  band  gap  was  increased  by  the  presence 

of  oxygen  vacancies.  From  this  argument  and  our  experimental  results,  it  appears  that 

our  crystal  was  treated  with  low  oxygen  pressure  (below  0.5atm)  and  therefore  the 

density  of  oxygen  vacancies  is  low  (see  Figure  2.12a).  The  largest  absorption 

coefficient  was  obtained  in  Figure  3.5.  In  this  configuration,  E  //  c  and  the  absorption 

coefficient  at  500nm  was  a  =2.32  cm'*  .  The  lowest  absorption  coefficient  was 

obtained  when  E  _L  c  in  Figure  3.4  and  absorption  coefficient  a  at  500nm  was  a  =0.88 

cm'*  .  The  band  gap  energy  can  be  found  from  the  absorption  spectrum.  The  band 

gap  energy  for  the  E  //c  case  is  E  =(hc/X  )  =  3.08eV  and  for  the  E  ±  c  case  it  is 

©  © 

3. 15eV. 

In  the  900-2500nm  wavelength  range  no  major  absorption  was  observed. 
Consequently,  the  absorption  spectrum  in  this  region  was  not  included.  As  seen  from 
Figure  3.7,  the  absorption  coefficient  was  zero  between  2.5ji  m  and  6.05jl  m,  and  there 
is  a  broad  absorption  band  between  8^1  m  and  15fi  m  where  the  absorption  is  quite 
high  (~  10  cm'*  ).  Last  [Ref.  37]  discussed  the  infrared  absorption  mechanism  in 
BaTiO-j  .  If  a  BaTiCA  crystal  contains  N  unit  cells,  and  each  unit  cell  has  five  atoms, 
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the  crystal  has  15N  degrees  of  freedom.  In  this  case,  3N  degrees  of  freedom  relates  to 
translational  motion  and  3N  to  torsional  motion  of  the  unit  cell.  The  remaining  9N 
modes  are  related  to  vibrational  degrees  of  freedom.  The  nine  vibration  modes  of  the 
unit  cell  can  be  classified  into  the  three  vibrations  of  Ba  against  the  TiO^  group  and  six 
internal  TiO^  vibrations.  TiO-j  group  was  considered  as  a  single  atom  situated  at  the  Ti 
position.  The  six  vibrations  of  Tio^  reduce  to  two  IR  active  vibrations. 


Figure  3.8  IR  active  normal  vibrations  of  a  TiO^  octahedron  [Ref.33]. 

In  Figure  3.8,  the  high  frequency  stretching  Vj  and  lower  frequency  bending  V2 
vibrations  are  illustrated.  The  stretching  vibration  is  due  to  a  change  in  length  of  Ti-Oj 
bond  and  the  bending  vibration  is  due  to  a  change  in  the  -Ti-Oj  bond  angle.  Then 
vibrational  motions  are  primarily  motions  of  oxygen  ions.  Last  [Ref.  37]  reported  two 
absorption  bands,  one  between  15-20^  m  and  another  one  around  20-30p  m  for  cubic 
BaTiOj  .  He  concluded  that  the  15-20^1  m  band  is  due  to  stretching  vibration  and 
20-30n  m  band  is  due  to  bending  vibration.  For  the  tetragonal  unit  cell,  he  observed 
only  the  stretching  vibration  at  9-20ji  m.  Our  results  showed  a  broad  absorption  band 
at  8-15n  m.  This  band  is  at  higher  frequencies  than  Last's.  This  may  be  due  to  the 
poling  process  performed  during  the  crystal  manufactured.  Poling  decreased  the  Ti-Oj 
distance  which  increased  the  vibrational  frequency. 

B.  PHOTOCURRENT  MEASUREMENT  EXPERIMENTS 
1.  Purpose  of  this  experiment 

Photocurrent  measurements  were  made  to  determine  the  pyroelectric  effect 
and  the  bulk  photovoltaic  effect  (BPVE).  In  order  to  determine  the  quantum  efficiency 
in  BaTiO^  ,  we  needed  to  measure  the  pyroelectric  current  and  the  short  circuit 
photovoltaic  current  (SCPVC). 
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2.  Procedure 

During  these  experiments,  a  Lexel  Model  95  Ar+  laser  was  used  as  a  coherent 
light  source.  It  was  operating  on  a  single  transverse  mode  but  on  multi  longitudinal 
modes.  For  the  photocurrent  measurements,  a  single  longitudinal  mode  laser  light  was 
not  vital.  The  output  beam  from  the  laser  was  1.3mm  diameter  and  vertically  polarized. 
In  order  to  avoid  the  photon  drag  effect,  platinum  electrodes  (the  area  of  the  electrodes 
was  0.25cm"  )  were  attached  to  the  polar  faces  of  the  crystal  which  were  perpendicular 
to  the  c  axis  and  parallel  to  the  laser  beam  (see  Figure  3.9). 


The  crystal  was  mounted  on  a  protoboard  and  this  system  was  shielded  by  an 
aluminum  box  in  order  to  reduce  environmental  interference.  In  this  setup,  the  c  axis  of 
the  crystal  was  always  parallel  to  the  electric  field  of  the  incident  beam.  The  short 
circuit  current  measurements  were  made  with  a  Keithley  617  electrometer  and  its 
output  feed  a  Hawlett-Packard  7090A  plotter.  These  experiments  were  made  at  night 
time  in  the  dark,  in  order  to  avoid  stray  light  from  the  environment.  For  convenience, 
200sec  duration  light  pulses  were  used.  A  black  card  was  used  to  chop  the  incident 
laser  light.  Photo-induced  currents  were  measured  at  three  different  wavelengths 
(514.Snm,  4S8nm.  457.9nm)  and  three  different  incident  laser  power  (25mW,  50mW, 
75mW)  at  each  wavelength  for  two  different  configurations.  The  measurements  were 
made  at  a  room  temperature  of  23°C. 
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3.  Results  of  the  experiments 


Figure  3.10  Measured  photocurrent  along  +c  axis. 


59 


TIME  sec 


Figure  3.11  Measured  photocurrent  along  -c  axis. 
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Figure  3.12 


i  (the  pyroelectric  current)  versus  incident  power. 
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Figure  3.13  iph  (SCPVC)  versus  incident  power. 
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4.  Interpretation  of  the  results 

Immediately,  on  unblocking  the  laser  beam  (OX  positions  in  Figure  3.10  and 
3.11),  a  high  peak  photocurrent  ( ~  100-200pA)  was  observed.  The  current,  after  the 
high  peak  showed  a  slow  decay  (50sec),  eventually  reaching  a  steady- state  value  of 
l-9pA.  On  blocking  the  laser  beam  (OF  position  in  Figure  3.10  and  3.11),  there  was  an 
immediate  transient  current,  1.17  times  larger  than  the  initial  transient  current  in 
magnitude  and  direction  but  in  the  opposite  direction.  This  time  the  current  decayed 
slowly  to  zero.  We  also  observed  that,  when  the  c-axis  direction  was  reversed,  the 
polarity  of  the  photo-induced  currents  reversed  (see  Figure  3.11).  If  we  compare  Figure 
3.10  and  3.11,  we  see  these  differences:  1)  The  polarities  of  the  photocurrents  are 
different.  2)  The  amplitude  of  the  photocurrents  in  Figure  3.11  is  1.6  + /-0. 3  times 
larger  than  the  amplitude  of  the  photocurrents  in  Figure  3.10.  This  is  probably, 
because  of  absorption.  The  absorption  coefficient  in  Figure  3.5  is  1.3  times  larger  than 
the  absoption  the  absoprtion  coefficient  in  Figure  3.6.  This  comparison  suggests  that 
the  the  higher  absorption  coefficient  produced  the  higher  peak  current  (see  Figure  3.5 
and  3.11).  The  large  transient  currents  are  pyroelectric  in  origin  and  the  steady-state 
currents  are  BPVE  in  origin  as  will  be  explained  below. 

TRANSIENT  CURRENT:  When  the  laser  beam  is  unblocked,  the  crystal 
temperature  rises  (~  10°C),  initially  at  the  rate  {(dT/dt)>0).  This  causes  a  pyroelectric 
current  to  flow,  obeying  Equation  2.12,  along  the  c-axis  al  the  crystal.  The  rise  time  of 
this  current  was  ~  200fi  sec  (see  Figure  3.10  and  3.11).  This  can  not  be  seen  from  the 
figures  because  the  plotter  response  time  was  too  slow.  It  was  measured  by  a  digital 
Nicolet  oscilloscope.  The  rise-time  was  determined  by  the  crystal  capacitance  ( ~  25pF) 
and  the  crystal's  internal  RC  equivalent  circuit  [Ref.  10].  As  the  temperature  increased 
to  its  new  steady-state  value,  (dT/dt)  decreased.  This  caused  the  pyroelectric  current  to 
decay  with  a  thermal  relaxation  time  (40sec)  to  a  steady-state  value,  which  exists  as 
long  as  the  crystal  was  illuminated.  In  principle,  when  the  crystal  is  at  thermal 
equilibrium  at  the  new  temperature,  then  (dT/dt)  was  zero.  This  means  that  the 
pyroelectric  current  vanishes.  But  the  steady-state  current  is  not  always  zero  (see  next 
section).  When  the  laser  beam  was  blocked,  the  crystal  started  to  cool  down.  At  the 
beginning,  its  rate  of  cooling  was  maximum  and  then  decayed  slowly  to  zero. 
Eventually  the  crystal  came  into  thermal  equilibrium  with  the  room  temperature.  This 
cooling  down  process  is  reflected  in  the  current,  which  is  in  the  direction  opposite  to 
the  initial  warming,  because  ((dT/dt)  <0}. 
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As  mentioned  earlier,  photocurrent  measurements  were  made  at  three 
wavelengths  (514.5nm,  488nm,  457. 9nm)  and  three  different  incident  beam  powers  at 
each  wavelength  for  two  different  crystal  orientations.  Figure  3.6  shows  the  peak 
transient  current  versus  the  incident  beam  power.  We  see  that  the  transient  current  is 
linearly  dependent  on  the  incident  beam  power  and  it  has  wavelength  dependence.  In 
Figure  3.12,  slopes  of  the  lines  are:  0.8  pA/mW  at  X  =514.5nm,  1.2pA/mW  at  X 
=  4SSnm  and  1.6pA/m\V  at  X  =457.9nm. 

STEADY-STATE  CURRENT:  This  current  is  due  to  the  bulk  photovoltaic 
effect  (BPVE)  and  is  a  result  of  anisotropic  charge  transport  (see  Section  1I-C).  It 
existed  as  long  as  the  crystal  was  illuminated  by  the  laser  beam.  It  shows  up  in  Figure 
3.6  and  3.7  as  superimposed  on  the  transient  current  (see  Figure  3.10  and  3.11). 

As  we  see  from  Figure  3.13,  the  short  circuit  photovoltaic  current  (SCPVC)  is 
linearly  dependent  on  intensity.  Thus,  this  experimental  observation  verifies  Equation 
2.18  given  in  the  thoery  section.  In  Figure  3.13,  slopes  of  the  lines  are:  0.044pA/mW 
at  X  =  514.5nm,  0.06  at  X  =  4SSnm  and  0.11  at  X  =  457.9nm.  The  SCPVC  increased 
with  frequency,  which  is  consistent  with  the  absorption  coefficients  in  Equation  2.18. 
Previous  experiments  showed  that  the  absorption  coefficient  increases  with  frequency. 
The  SCPVC  is  only  observed  along  the  c-axis.  This  means  that  the  field  along  the  c 
axis  of  the  crystal  separate  carriers.  This  statement  is  in  complete  agreement  with 
Equation  2.17.  In  Figure  3.11,  when  the  c-axis  was  reversed,  the  polarity  of  the 
SCPVC  reversed  (see  Equation  2.15). 
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C.  BEAM-COUPLING  EXPERIMENTS 

1.  Purpose  of  this  experiment 

In  order  to  calculate  the  quantum  efficiency  of  BaTiO^  ,  we  needed  to  know 
the  charge  transport  length.  We  know  that,  theoretically,  charge  migration  between 
traps  cause,  approximately  a  90°  phase  shift  between  the  hologram  grating  (index  of 
refraction  distribution)  and  the  interference  pattern  (see  Section  11-G).  This  phase 
angle  and  the  charge  transport  length  can  be  calculated  using  the  results  of  following 
experimental  procedure. 

2.  Procedure 

A  Lexel  Model  95  Ar+  laser  was  used  as  a  coherent  light  source.  The  output 
beam  from  this  laser  was  vertically  polarized.  We  needed  a  horizantally  polarized  beam 
for  our  purposes.  Two  quarter  wave  plates  were  used  to  accomplish  this.  As  we  see 
from  Figure  3.14,  beam  splitter  (BS)  split  the  output  beam  into  two  equal  intensity 
beams.  One  beam  followed  the  path  pi  (36cm)  and  mirror  m2  reflected  it  back  to  the 
crystal,  at  an  angle  9  =  17.5°  with  respect  to  the  crystal  surface  normal.  Mirror  ml, 
reflected  the  second  beam  onto  mirror  m3.  This  mirror  reflected  the  second  beam  onto 
the  crystal,  making  angle  9  =  17.5°  with  the  surface  normal.  The  second  beam  path  p2 
was  33cm  long.  Two  beams  crossed  each  other  at  an  angle  ~  7°  in  the  crystal  (recall 
Snell's  law).  It  wTas  important  to  keep  the  difference  in  the  path  lengths  (in  our  case 
pl-p2  =  3cm)  of  the  twTo  beams  within  the  coherence  length  limit  of  the  laser,  which 
was  about  ~  6cm.  Otherwise  the  interference  pattern  inside  the  crystal  could  not 
occur.  The  intensity  of  the  first  beam  was  1/2  of  the  second  beam  (A  =  0.5).  Initially, 
the  two  beams  were  kept  on  the  crystal  until  the  grating  was  established.  Then  the 
strong  beam  wras  blocked  to  determine  the  energy  transfer  from  the  strong  beam  to  the 
weak  beam. 
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Figure  3.14  Experimental  set-up  for  beam-coupling.. 
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3.  Results  of  the  experiments 


Figure  3.15  The  fractional  energy  transfer  (F)  at  X  =488nm. 

The  numbers  on  the  curves  denote  the  fractional  energy  transfer  (F).  The  I  is  the 
intensity  of  the  signal  beam  when  the  reference  beam  was  turned  off.  The  I  ^  is  the 
intensity  of  the  signal  beam  when  the  reference  beam  was  turned  on. 
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Figure  3.16  The  fractional  energy  transfer  (F)  at  X  =  457.9nm. 

The  numbers  on  the  curves  denote  the  fractional  energy  transfer  (F).  The  I  ^  is  the 
intensity  of  the  signal  beam  when  the  reference  beam  was  turned  on.  The  I  is  the 
intensity  of  the  signal  beam  when  the  reference  beam  was  turned  off. 

4.  Interpretation  of  the  results 

As  is  shown  in  Figure  3.15  and  3.16,  when  the  strong  beam  was  blocked  (OFF 
conditions  in  the  figures),  the  intensity  of  the  weak  beam  decreased  and  reached  a 
steady-state  value  I  .  After  10  seconds,  the  strong  beam  was  unblocked  (ON  positions 
in  the  figures),  and  the  intensity  of  the  weak  beam  increased  and  reached  a  steady-state 
value  I  ^  .  This  is  the  manifestation  of  energy  transfer  in  beam-coupling.  Energy 
transfer  occurs,  because  the  refractive  index  grating  (hologram  grating)  scatters  the  two 
beams  that  produce  it.  Using  Equation  2.53  and  the  data  obtained  from  the  beam¬ 
coupling  experiments,  the  diffraction  efficiency  of  the  hologram  grating  was  calculated. 
Figure  3.17  shows  the  diffraction  efficiency  versus  wavelength  of  the  laser  beam.  The 
direction  of  energy  transfer  is  determined  by  the  direction  of  the  c-axis  of  the  crystal 
and  the  phase  shift  angle.  According  to  Valley  [Ref.  24],  the  direction  of  energy  transfer 
depends  on  the  type  of  carrier  which  dominates  in  charge  transport.  In  our 
experimental  set-up,  the  strong  beam  enters  the  crystal  at  the  -c-axis  edge  and  beam- 
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Figure  3.17  Diffraction  efficiency  vs  wavelength. 
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coupling  occurs  in  the  +c  direction.  This  means  that  in  our  crystal,  the  majority 
carriers  are  electrons  (according  to  Valley's  theory’).  Our  observations  agree  with 
Gunter's  [Ref.  3S]  arguments.  These  arguments  are:  l)  Energy  transfer  occurs  only  if 
the  two  writing  beams  have  different  intensities.  In  our  experiment,  A=  Is  / If  =0.5.  2) 
Energy  is  always  transferred  from  the  strong  beam  to  the  weak  beam.  3)  If  |(pQ  -<p  |  >  tt 
2  ,  where  (pQ  is  the  phase  difference  between  the  two  beams  and  (p  is  the  phase  shift  of 
the  index  grating  from  the  intensity  pattern,  then  the  energy  transfer  is  oscillatory.  In 
our  experiments,  the  difference  between  the  path  lengths  of  the  two  beams  was  3cm. 
This  distance  introduced  a  phase  difference  between  the  two  incident  beams.  This  is 
why,  we  observed  an  oscillatory’  response  under  steady-state  conditions  in  Figure  3.15 
and  3.16.  As  the  path  difference  increased,  the  oscillations  increased. 

D.  CALCULATION  OF  THE  QUANTUM  EFFICIENCY  OF  BATIO3 

In  this  section,  the  quantum  efficiency  of  BaTiO^  will  be  calculated,  using  the 
experimental  data  and  the  equations  derived  in  Section  II-H.  In  Section  III-B,  the 
pyroelectric  current  was  given,  this  data  will  be  used  to  calculate  the  photorefractive 
absorption  coefficient.  Using  Equation  2.66  and  the  measured  value  of  the  pyroelectric 
current  we  found  the  photorefractive  absorption  coefficient  ap  .  For  example,  at  X 
=  4SSnm,  the  measured  value  of  the  pyroelectric  current  was  ipy  =  27.26pA  (see  Figure 
3.12).  Using  Equation  2.66,  the  photorefractive  absorption  coefficient  was  found  to  be 
ap  =0.00S6cm'^  for  the  beam  power  25m\V  at  X  =  4SSnm.  Similar  calculations  were 
made  at  X  =514.5nm  and  X  =457.9nm  for  three  different  incident  beam  powers 
(25m\V,  50mW,  75mW).  Figure  3.17  shows  the  photorefractive  absorption  coefficient 
versus  wavelength  of  the  laser  beam.  The  absorption  coefficients  were  measured  by  the 
spectrophotometer  (for  example,  see  Figure  3.3)  have  the  same  wavelength  dependence 
with  Figure  3.17.  The  comparison  between  the  a  and  the  ap  shows  that  the 
photorefractive  absorption  coefficient  is  approximately  200  times  smaller  than  the 
absorption  coefficients  were  found  in  Section  III-A.  For  example,  from  the 
spectrophometric  measurement  the  absorption  coefficient  was  found  to  be  «  =2cirf^ 
at  X  =  457.9nm.  The  photorefractive  absorption  coefficient  from  the  pyroelectric 
current  measurements  wras  found  to  be  ap  =0.013cm'^  at  X  =457.9nm.  Moreover, 
since  Glass  [Ref  39]  has  shown  that  the  photo-induced  change  in  the  refractive  index  is 
strongly  dependent  on  the  incident  beam  intensity  and  since  ap  was  measured  here  is 
dependent  on  the  incident  beam  power,  w’e  have  clear  causal  connection  between  our 
measured  ap  and  photorefraction. 
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As  stated  in  Section  II-B,  the  anisotropy  constant  K  is  the  measure  of  the 
difference  between  the  probabilities  of  exicitation  and  recombination  along  the  +  c  and 
-c  axes  of  the  crystal.  The  anisotropy  constant  is  needed  to  calculate  the  quantum 
efficiency.  Using  Equation  2.68  and  the  experimentally  measured  values  of  the  SCPVC 
(see  Figure  3.13),  the  anisotropy  constant  were  calculated.  For  example,  at  X  =488nm 
the  measured  value  of  the  absorbed  power  was  O  =  2. 15*  lO-*^  Watt/cm,  using  i 
=  27.26pA  in  Equation  2.66  and  a  measured  value  of  SCPVC  was  iph  =  1.62pA  (see 
Figure  3.13).  Putting  these  numbers  in  Equation  2.68  the  anisotropy  constant  was 
found  to  be  K  =0.39*10*9  A  cm  /Watt.  Similar  calculations  were  made  at  X 
=  514.5nm  and  X  =457.9nm  for  three  different  incident  beam  powers  (25mW,  50mW, 
75mW).  In  Figure  3.18,  the  anisotropy  constant  is  plotted  versus  wavelength.  The 
phase  shift  angle  (pph  is  the  last  parameter  involved  in  the  calculation  of  the  quantum 
efficiency.  Using  the  measured  value  of  the  fractional  energy  transfer  (see  Figure  3.15 
and  3.16),  M  (from  Equation  2.49)  and  A  =  0.5  in  Equation  2.70,  <pph  was  calculated. 
At  X  =488nm,  the  measured  value  of  the  fractional  energy  transfer  was  F  =  0.22  (see 
Figure  3.15),  M=  12.57  (from  Equation  2.49)  and  the  intensity  ratio  of  the  signal  and 
the  reference  beams  was  A  =  0.5.  Using  Equation  2.70,  the  phase  shift  angle  was  found 
to  be  (pph  =52.1°.  Similar  calculations  were  made  at  X  =514.5nm  and  at  X 
=  457.9nm.  Figure  3.20  shows  the  phase  shift  angle  (pph  versus  wavelength.  Using  the 
values  of  the  anisotropy  constant  (from  Figure  3.19),  the  phase  shift  angle  (from  Figure 
3.20)  and  the  wave  number  of  the  the  hologram  grating  (K  =  4*;t  sinG  /X  )  in  Equation 
2.72,  the  quantum  efficiency  was  calculated.  For  example,  at  X  =488nm,  the  value  of 
the  anisotropy  constant  wras  K  =0.39*10*9  ^  cm/W,  the  value  of  the  phase  shift  angle 
wras  52.1°  and  the  wrave  number  of  the  hologram  grating  wras  K  =  8.25*10^  m.  Using 
Equation  2.72,  the  quantum  efficiency  of  BaTiO^  wras  found  to  be  5.96*10*^  .  Similar 
calculations  were  made  at  X  =514.5nm  and  at  X  =457.9nm.  Figure  3.21  plots  the 
quantum  efficiency  is  plotted  versus  wavelength. 
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Figure  3.18  Photorefractive  absorption  coefficient  of  BaTiOj. 
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Figure  3.19  Anisotropy  constant  of  BaTiCA. 


73 


WAVELENGTH  (NM) 


Figure  3.20  Phase  shift  angle  between  gratings. 
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QUANTUM  EFFICIENCY 


Figure  3.21  Quantum  efFiciency  of  BaTiO-j. 
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E.  DISCUSSIONS 


This  section  summarizes  the  photorefractive  effect  in  BaTiO^  quantitatively. 
This  will  be  done  with  a  two  beam-coupling  example  at  one  wavelength.  At  X 
=  457.9nm,  the  wave  number  of  the  intensity  interference  pattern  was  K=  8.25*  10^  m 
(K  =  4*rc  sinG  X  ,  where  0  =  17.5°  is  angle  with  crystal  surface  normal).  If  3.9*10*^ 
photons  per  sec  were  absorbed  (from  Equation  2.63  for  the  value  of  the  absorbed 
power  0  =  1.72*10*^  Watt  cm),  photo-carriers  (in  our  case  electrons  are  majority)  will 
be  ,.eed  from  the  absorption  centers  (see  Section  II-E)  Under  the  influence  of  BPVE 
(see  Section  II-H),  carriers  will  migrate  from  the  illuminated  regions  to  the  dark 
regions.  This  process  sets-up  a  space  charge  field,  since  we  assumed  a  one  dimensional 
charge  transport  (c-axis  direction)  in  Section  II-C.  Using  the  measured  value  of  the 
fractional  energy  transfer  F  =  0. 17  (see  Figure  3.15),  A  =  0.5  (the  beam  ratios)  and 
M  =  10.06  (from  Equation  2.49)  in  Equation  2.70,  the  phase  shift  angle  was  found  to  be 
<p  h  =  23.72°.  With  this  value  of  the  phase  shift  angle  in  Equation  2.71,  the  charge 
transport  length  was  found  to  be  53nm.  Using  the  measured  value  of  the  fractional 
energy  transfer  F  =  0.17  in  Equation  2.53,  the  diffraction  efficiency  of  the  hologram 
grating  was  found  to  be  q  =0.0072.  Using  this  value  of  diffraction  efficiency  and  a 
3.6mm  hologram  grating  (l=t,cos0  )  in  Equation  2.49,  the  amplitude  of  refractive 
index  change  was  found  to  be  5  n^  =2.45*10'^  .  Using  from  Table  1  the  values  of 
the  extraordinary  index  of  refraction  n^  =  2.365  and  the  electrooptic  coefficient  along 
the  c-axis  of  the  crystal  r^  =23*10  m/V  and  the  value  of  the  refractive  index 
change  in  Equation  2.41,  the  magnitude  of  the  space  charge  field  was  found  to  be  16.11 
kV.  m.  The  quantum  efficiency  of  BaTiO^  at  this  wavelength  was  found  to  be  2.13*10""* 
(see  Figure  3.21).  Using  this  value  of  the  quantum  efficiency  in  Equation  2.29,  the  the 
generation  rate  of  electrons  was  8.3*  10^  per  sec.  If  we  consider  the  density  of 
photorefractive  center  as  10*^  cm'-*  [Ref.  19],  then  in  a  volume  of  6.6*10'-*  cm'*  (Ab 
=  0.0132cm“  ancFt  =  0.5cm)  6.6*10*-*  electrons  were  found.  Glass  et  al.  [Ref.  31]  found 
the  quantum  efficiency  of  LiNbCU  as  0.17,  this  value  of  the  quantum  efficiency  gave  a 
electron  density  that  was  three  times  the  density  of  the  photorefractive  center,  which  is 
impossible.  Our  calculation  of  the  quantum  efficiency  gave  a  value  of  the  electron 
density  (see  above  photo-excited  electron  density  is  smaller  than  the  impurity  density) 
that  is  consistent  with  the  density  of  oxygen  vacancies. 
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IV.  CONCLUSIONS 


The  fact  that  the  absorption  coefficients,  measured  in  the  visible  region  of  the 
spectrum,  were  found  to  be  different  for  each  of  the  six  crystal  orientations  shows  that 
the  absorption  centers  are  distributed  anisotropically  in  the  crystal.  From  these 
measurements,  the  band  gap  of  BaTiCU  was  found  to  be  E  (1)  =3.08eV  and  E  (2)  = 
3.l5eV.  A  broad  plateau  was  observed  between  880nm  and  the  band  edge  at  460nm. 
This  was  due  to  the  low  oxygen  pressure  used  during  crystal  treatment.  In  the  infrared 
region  of  the  spectrum,  a  broad  band  8-154  m  was  observed  and  this  comes  from 
stretching  vibrations  of  TiC>3  atoms.  At  X  =500nm  the  maximum  absorption 
coefficient  was  a  =2. 32cm' ^  and  the  minimum  absorption  coefficient  was  a 
=  0.88cm'l  .  These  numbers  are  too  large  to  be  accounted  by  the  photorefractive 
absorption  centers  alone.  Since  Glass  [Ref.  39|  has  shown  that  the  photo-induced 
change  in  the  refraction  index  is  strongly  dependent  on  the  incident  beam  intensity  and 
since  our  measurements  of  the  photorefractive  absorption  coefficient  ap  were 
dependent  on  the  incident  beam  power,  this  suggests  the  causal  connection  between 
our  measuement  of  ap  and  photorefraction  (see  Figure  4.1  and  3.18). 


Figure  4.1  Intensity  vs  the  index  of  refraction  change  [Ref.39]. 

The  photorefractive  absorption  coefficient  were  determined  from  the  pyroelectric 
current  measurements.  The  maximum  value  of  the  photorefractive  absorption 
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coefficient  was  found  to  be  ap  =  0.0145cm'*  at  X  =457.9nm  for  a  beam  intensity  of 
75mW.  The  minimum  value  of  a  was  found  to  be  0.0086cm'*  at  X  =488nm  for  a 

P 

beam  intensity  of  25mW.  We  see  that  the  photorefractive  absorption  coefficient  is  ~ 

200  times  smaller  than  the  absorption  coefficient  obtained  from  the  spectrophotometric 

measurements.  This  important  point  has  been  overlooked  previously  in  the  literature. 

From  the  short  circuit  photovoltaic  current  (SCPVC)  and  the  absorption  coefficients 

derived  from  the  pyroelectric  current  measurements,  the  anisotropy  constant  was  found 

to  be  K  =  0.4*  1  O''*  (Ampere. cm/watt)  at  X  =457.9nm.  With  this  information  the 

quantum  efficiency  times  charge  transport  length  was  found  to  be  QLph  =  1.13*10'^ 

cm.  From  beam-coupling  experiments,  the  amount  of  fractional  energy  transfer 

between  two  beams  was  found  to  be  ~  F  =  0.18.  Combining  this  information  with  the 

results  of  coupled  wave  theory,  the  diffraction  efficiency  of  a  hologram  grating  was 

found  to  be  n  =0.018  at  X  =514.5nm.  The  maximum  value  of  the  phase  shift  angle 

between  the  hologram  grating  and  the  intensity  pattern  at  X  =514.5nm  was  78.9°.  The 

minimum  charge  transport  length  in  photovoltaic  process  was  Lph  =  53nm.  The 

quantum  efficiency  increases  with  the  frequency  of  incident  light.  The  maximum  value 

1 

of  the  quantum  efficiency  was  Q=  2.13*10  (at  X  =457.9nm),  and  the  minimum  value 
of  the  quantum  efficiency  was  Q=  1.1 3*  10°  (at  X  =514.5nm).  For  absorption  of 
3.9*10*“*  photons  per  sec  at  X  =457.9nm  and  with  this  value  of  the  quantum 
efficiency,  the  generation  rate  of  electrons  was  8.3*10***  per  sec.  If  the  density  of 
absorption  centres  is  10***  cm'^  ,  then  in  the  volume  occupied  by  the  beam  in  the 
crystal  (6.6*10'"*  cm"*  ),  there  are  6.6*10*"*  electrons.  Young  et  al.  [Ref.  32]  found  the 
quantum  efficiency  of  LiNbO"*  as  Q=2*10'"*  and  the  charge  transport  length  as  Lph 
=  50tim  (wavelength  was  not  given  in  the  literature).  But  they  used  the  large  values  for 
the  absorption  coefficient.  This  compares  well  the  quantum  efficiency  of  BaTiO^  as 
Q=  2.13*10'“*  and  the  charge  transport  length  as  Lph  =53nm  at  X  =457. 9nm  that  we 
found.  Since  the  lifetimes  of  the  carriers  are  still  unknown,  the  Dember  effect  can  be 
used  to  determine  the  lifetimes  in  future  measurements. 


78 


LIST  OF  REFERENCES 


1.  Chen,  F.S.,  J.AppLPhys.,  40,  3389,  ( 1969) 

2.  Kukhtarev,  N.V.,  Markov,  V.B.,  Odulov,  S.G.,  Soskin,  M.S.,  Vinetskii,  V.L 
Ferroelectrics,  22,  949,  ( 197^) 

3.  Jona,  F.,  Shirane,  G.,  Ferroelectric  Crystals,  Pergamon  Press,  (1962) 

4.  Kittel,  C.,  Introduction  to  Solid  State  Physics,  Sixth  Edition,  Wiley  (1986) 

5.  Devonshire,  A.F.,  Phil.Mag.,  40,  1040,  (1949) 

6.  Slater,  J.C.,  Phys.  Rev.,  78,  748,  (1950) 

7.  Merz,  W.J.,  Phys.  Rev.,  76,  1221,  (1949) 

8.  Chynoweth,  A.G.,  J.Appl.Phys.,  27,  78,  (1956) 

9.  Eustace,  L.D.,  Devon,  G.C.,  Optical  Radiation  Detectors,  Wiley  (1984) 

10.  Chynoweth,  A.G.,  Phys. Rev.,  102,  705,  (1956) 

11.  Glass,  A.  M.,  Appl.Phys.Lett.,  25,  233,  (1974) 

12.  Gower,  M.C.,  Opt.Lett.,  11,  458,  (1986) 

13.  Brody,  P.  S.,  Sol. St. Comm.,  12,  673,  (1973) 

14.  Koch.  W.T.H.,  Munser,  R.,  Ruppel,  W.,  Wurfel,  W.  Sol. St. Comm.,  17,  847, 

\iy  I j) 

15.  Bube,  R.H.,  Photoconductivity  of  Solids,  Wiley  (1960) 

16.  Glass,  A.M.,  Chanussot,  G.,  Phys. Lett  ,  59A,  405,  (1976) 

17.  Baltz,  R.Von.,  Kraut,  W.,  Sol.St.Comm.,  26,  961,  (1978) 

18.  Klein,  M.B.,  Schwartz,  R.N.,  J.Opt.Soc.AmB.  ,  3,  293,  (1986) 

19.  Feinberg,  J.,  Ducharme,  S.,  J.Opt.Soc.AmB.  3,  283,  (1986) 


79 


20. 

21. 


Hellwarth,  R.V.,  Strohkendl,  F.P.,  Jonathan,  J.M.C.,  Opt.Lett.,  15,  313,  (1986) 

Young,  L.,  Moharam,  M.G.,  J.Appl.Phys.,  48,  3230,  (1977) 

Feinberg.  J.,  Heiman,  D.,  Tanguav,  A.R.,  Helwarth,  RAV.  J.Appl.Phys.,  51, 
1297, ( 1980)  ff  J  -  - 


23.  Young,  L.,  Moharam,  M.G.M.,  Gaylord,  T.K.,  Magnusson,  R.  Ferroelectrics, 

27,253,(1^80)  6 

24.  Valley,  G.C.,  J.Appl.Phys.,  59,  3363,  (1986) 

25.  Yariv,  A.,  Quantum  Electronics,  Second  Edition,  Wiley,  (1985) 

26.  Johnston,  D.,  J.Appl.Phys.,  41,  3279,  (1970) 

27.  Vahey,  W.D.,  J.Appl.Phys.,  46,  3510,  (1975) 

28.  Young,  L.,  Moharam,  M.G.,  Gavlord,  T.K.,  Magnusson,  R.  J.Appl.Phvs.  50, 

5642.7197$)  ff 

29.  Kogelnik,  H.,  Bell.Syst.Techn.J.,  48,  2909,  (1969) 

30.  Amodei,  J.J.,  Staebler,  D.L.,  J.Appl.Phys.,  43,  1042,  (1972) 

31.  Glass,  A.M.,  Peterson,  G.E.,  Negran,  J.T.,  An  Introduction  to  Electrooptic 
Devices.  Kaminow,  I.P.,  Academic  Press.  325,  (1974) 

32.  Young,  L.,  Moharam,  M.G.,  Guibaly,  F.E.,  Lun,  E.,  J.Appl.Phys.  50,  4201, 
(1979J 

33.  Kratzig,  E.,  Orlowski,  R.,  Opt.Qm.Eln.,  12,  495,  (1980) 

34.  Young,  L.,  Moharam,  M.G.,  Guibaly,  F.E.,  Lun,  E.  J.Appl.Phys.,  50,  4201, 


35.  Vahey,  W.D.,  J.Appl.Phys.,  46,  3510,  (1975) 

36.  Casella  R.C.,  Keller,  S.P.,  Phys.Rev.,  15,  1469,  (1959) 

37.  Last,  J.T.,  Phys.Rev.,  105,  1740,  (1956) 

38.  Gunter,  P.,  Phys.Rep.,  93,  199,  (1982) 

39.  Glass,  A.M.,  Opt.Eng.,  17,  470,  (1980) 


80 


INITIAL  DISTRIBUTION  LIST 


No.  Copies 

1.  Defense  Technical  Information  Center  2 

Cameron  Station 

Alexandria,  Virginia  22304-6145 

2.  Library.  Code  0142  2 

Naval  "Postgraduate  School 

Monterey,  California  93943-5000 

3.  Donald  L.  Walters  10 

Code  61  We 

Naval  Postgraduate  School 
Monterey,  California  93943 

4.  Nathaniel  E.  Glass  ^ 

Code  61  Gc 

Naval  Postgraduate  School 
Monterey,  California  93943 

5.  Gazanfer  Unal  4 

Kaptan  Arif  Cad  Coban  Yildizi  Sok 

No  =  26/19  Erenkoy 
Istanbul  Turkey 

6.  Deniz  Harp  Okulu  Kutuphanesi  1 

Tuzla 

Istanbul  Turkey 

7.  Bogazici  Universitesi  1 

Kutuphane 

Bebek 

Istanbul  Turkey 


81 


(TV 


3 


f 


1 8070  2 


-•  l 


Thesis 
U353 
c.  1 


OuOS 


Unal 

Measurements  of  the 
microscopic  properties 
contributing  to  photore 
fraction  in  barium  tita 
nate . 


thesU353 

Measurements  of  the  microscopic  properti 


i'ii  mu  n  iiiiii  iiiiiiiiiiiiiu  ii  ii  i  iiiii  mu  mu  i  ii 

3  2768  000  70972  9 
DUDLEY  KNOX  LIBRARY 


